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Abstract 

We develop the theory of Abelian functions associated with algebraic curves. The growth in computer 
power and an advancement of efficient symbolic computation techniques has allowed for recent progress 
in this area. In this paper we focus on the genus three cases, comparing the two canonical classes of 
hyperelliptic and trigonal curves. We present new addition formulae, derive bases for the spaces of 
Abelian functions and discuss the differential equations such functions satisfy. 

1 Introduction 



In this paper we present addition formulae and differential equations satisfied by the Abelian functions with 
poles along the standard theta divisor, associated with the hyperelliptic and trigonal curves of genus three, 
i-^ ' Alongside the presentation of the new results, this paper offers a comparison of the formulae for the two 

C$ . canonical types of genus three curves. 

The most important new results are the addition formulae presented in Theorems 15.31 and 15.41 The first 
of these was obtained by following the ideas in [13] while the second requires the explicit derivation of a basis 
for the vector space of Abelian functions having poles of order at most four along the standard theta divisor. 
They both follow from the recent work in |15| which introduced new addition formulae for the Weierstrass p 
' and (T-functions along with generalisations to genus two. We also present and discuss the sets of differential 

equations satisfied by the functions. These results can be seen as a continuation of the work started in [T3] 
CN ■ and [1] using new efficient computational techniques first introduced in [18j . 

Curves of genus three can be categorised as either hyperelliptic or trigonal. We study two canonical 
examples, the (2,7) and the (3,4)-curves which are hyperelliptic and trigonal respectively. We define these 
terms and the curves formally in the next section. We develop the results of [T3] for the (3,4)-curve, 
completing the key sets of differential equations and deriving new addition formulae. We compare these 
results with the corresponding genus three hyperelliptic results for the (2, 7)-curve, furthering the results in 
[8] for this case. 

Many of the results we present can be viewed as generalisations of classic results for elliptic functions. 
' We will conclude the introduction below by reminding the reader of the relevant results from the theory 

of Weierstrass functions. Then in Section [2] we give the definitions of the Abelian functions, discussing the 
general properties they satisfy. We proceed in Section [3] to consider the problem of determining bases for 
the vector spaces of such functions, presenting explicit constructions of these for the genus three curves. In 
Section [5] we derive sets of differential equations satisfied by the functions and compare the two canonical 
genus three cases. Finally is Section [5] we present the new addition formulae. 

The classical results for Weierstrass elliptic function form a template for our theory. Let p(u) be the 
Weierstrass p-function, which as an elliptic function has two complex periods wi,W2' 

p(u + uj\) = p(u + uj 2 ) — p(u), for all u e C (1) 

The p-function has the simplest possible pole structure for an elliptic function and satisfies a number of 
interesting properties. For example, it can be used to parametrise an elliptic curve, 

y 2 = 4x 3 -g 2 x-g 3 , (2) 

where g 2 and g 3 are constants. It also satisfies the following well-known differential equations, 

(p'(u)) 2 = 4p(u) 3 - g 2 p(u) - g 3 , (3) 
p"(u)=6p(u) 2 ~y 2 . (4) 
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Weierstrass introduced an auxiliary function, a(u), in his theory which satisfies 



d 2 

pO) = lo s [ (T ( U )] • ( 5 ) 

The cr-function plays a crucial role in the generalisation and in applications of the theory. It satisfies the 
following two term addition formula, 

a(u + v)a(u - v) 

a(u) 2 a{v) 2 =^~^- <® 

In this document we present generalisations of equations ([1]) — © for higher genus functions. 

Elliptic functions have been the subject of much study since their discovery and have been extensively 
used to enumerate solutions of non-linear wave equations. Recent times have seen a revival of interest in 
the theory of Abclian functions, which have multiple independent periods, and so generalise the elliptic 
functions. The periodicity property is usually defined in association with an underlying algebraic curve. 
These functions have been shown to solve differential equations arising in mathematical physics and have 
been used in a variety of applications. 



2 Constructing Abelian functions 

In this paper we study in detail the case of functions associated with genus three curves. However, much of 
the theory is applicable to a wider space of curves and so we will include some general definitions. 

Definition 2.1. For two coprime integers (n,s) with s > n we define an (n, s)-curve as a non-singular 
algebraic curve defined by f(x,y) = 0, where 

f(x,y)=y n +p 1 (x)y n - 1 + p 2 (a:)i/"~ 2 + • • • + Pn-i(x)y - p n (x). (7) 

Here x, y are complex variables and Pj(x) are polynomials in x of degree (at most) \Js/n\. We define a 
simple subclass of the curves by setting Pj{x) = for < j < n — 1. Such curves are then defined by 

f(x, y) = y n - (x s + Xs^x*- 1 + ■ ■ ■ + X x x + A ) (8) 

and are called cyclic (n, s)- curves. We follow tradition and denote the curve constants by Xj for the 
cyclic (n, s)-curves and by Hj for the general (n, s)-curves. Note that in the literature the word "cyclic" is 
sometimes replaced by "strictly" or "purely (n-gonal)" . Also, some authors use a different normalisation 
with the coefficient of X s chosen to be 4 instead ofl. It is simple to move between the different normalisations. 

We denote the surface defined by such a curve as C . The genus of C is given by 

0=!(n-l)(*-l) 

and the associated functions to be defined shortly will be multivariate with g variables, u = (ui, . . . ,u g ). 
For example, the elliptic curve in equation ^ is a (2,3)-curve and the Weierstrass a and p-functions depend 
upon a single variable u. 

Given this equation for the genus we see that it can only be three when (n, s) = (2,7) or (3,4). Hence 
the two curves we will investigate in this paper are as follows: 

y 2 + (/iix 3 + fi 3 x 2 + fisx + fi?)y — x 7 + [i 2 x 6 + ^x 5 + fi6X 4 + Ms^ 3 + tiiox 2 + Hux + /U14, (9) 
y 3 + (fiix + HA)y 2 + (fJ-2X 2 + tmx + (i$)y = x 4 + fi 3 x 3 + j_i e x 2 + ^ 9 x + fj,i 2 . (10) 

Their cyclic restrictions are given by: 

y 2 = x 7 + X 6 x 6 + X 5 x 5 + A 4 x 4 + X 3 x 3 + X 2 x 2 + Xyx + A , (11) 
y 3 = x 4 + X 3 x 3 + X 2 x 2 + X lX + X . (12) 

We define a weight for the theory, denoted by wt and called the Sato weight. We start with 

wt(x) = —n, wt(y) = — s (13) 
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and then choose the weights of the curve parameters to be such that the curve equation is homogeneous. 
We see that for cyclic curves this imposes wt(Aj) = — n(s — j) while for the non-cyclic curves © and (fTUj) 
the constants nj have been labeled with the absolute value of their weights. The weights of all other objects 
may then be derived uniquely and will be commented on throughout the paper. Note that all the objects in 
this paper have a definite Sato weight and all the equalities are homogeneous with respect to this weight. A 
more detailed discussion of the weight properties may be found for example in [18) . 

The (n, s)-curves with n = 2 are generally defined to be hyperelliptic curves. Klein developed an approach 
to generalise the Weierstrass p-function to Abelian functions associated with hyperelliptic curves as described 
in Baker's classic texts [3] and [S]. This approach has motivated the general definitions in [S] and [33] of 
what we now call Kleinian p-functions. It is the properties of these and the generalised cr- function that are 
our objects of study. 

In the last few years a good deal of progress has been made on the theory of Abelian functions associated 
to those (n, s) curve with n = 3, which we label trigonal curves. In [10j . the authors of [S] furthered 
their methods to the trigonal cases, obtaining realisations of the Jacobian variety and some key differential 
equations between the functions. More recently the two canonical trigonal cases of the (3,4) and (3,5)-curves 
have been examined in |13j and [5] respectively. Both papers explicitly construct the differentials on the 
curve, solve the Jacobi inversion problem and obtain differential equations between the p-functions. Some 
of the properties of higher genus trigonal curves have been explored in |16| . The class of (n, s)-curves with 
?i = 4 are defined as tetragonal curves and have been recently considered for the first time in |18) and [19] . 
The lowest genus tetragonal curve has genus six. 

We now discuss how to construct these functions, starting by choosing bases for the space of differential 
forms of the first kind. A standard basis may be constructed from the Weierstrass gap sequence associated 
to the curve (see for example [7]). For any (2, 7)-curve the basis is 

^ / dx xdx x 2 dx \ 
\fy(x,y)' f y (x,y)' fyix.y)) 1 

while for any (3, 4)-curve it is 

d f dx xdx ydx \ 
^fy{x,vY fy(x,y)' f y {x,y)J' 
We choose a symplectic basis in Hi(C, Z) of cycles (closed paths) upon the compact Riemann surface defined 
by C. We denote these by {ai, «2, 013, Pi, ft, ft} and ensure they have intersection numbers 

atj • aj =0, ft • ft = 0, on ■ ft = 

We introduce dr as a basis of meromorphic differentials which have their only pole at 00. We will not 
explicitly use these in this paper and so for an explicit construction we refer the reader to Section 2.2.2 in 
[8] for the (2, 7)-case and to equation (2.22) in [13] for the (3, 4)-case. These bases are derived alongside 
a fundamental differential of the second kind which plays an important role in the theory of p-functions. 
Again, we will not explicitly use the differential in this paper and so refer the reader to [7] for an example of 
its construction and role. 

We can now define the standard period matrices associated to the curve as follows. 

^' = (^)v=l,2,3 ^ = M-1A8 

We define the period lattice A formed from lo', ui" by 

A = {moj' + nuj" , m,n€Z 3 }. 
The functions we treat in this paper are defined upon C 3 with coordinates usually given as 

u = (ui, u 2 , u 3 ). (16) 

Note that any point it £ C 3 can be given by 




E 

i=l 



du. 
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where the Pi are points upon C. The period lattice A is a lattice in the space C 3 . Then the Jacobian variety 
of C is presented by C 3 /A, and is denoted by J. We define k as the modulo A map, 

K : C 3 ^J. (17) 

For k = 1, 2, . . . define 2lfe, the Abel map from the k-th symmetric product Sym fe (C) to J. 

2l fc : Sym fe (C) — » J 

(P 1 ,...,P k )^ ( [ 1 du+---+ f " du) (mod A), (18) 

where the P 2 ; are again points upon C. Denote the image of the fc-th Abel map by and define the k-th 
standard theta subset (often referred to as the fc-th strata) by 

Q[k] = w [k] y 

where [—1] means that 

[— 1](U1,U2,U3) = (— «1, —U2, ~U 3 ). 

When k = 1 the Abel map gives an embedding of the curve C upon which we define £ = a; - ™ as the local 
parameter at the origin, 2ti(oo). We can then express y and the basis of differentials using £ and integrate 
to give series expansions for u. We can check the weights of u from these and see that they are prescribed 
by the Weierstrass gap sequence. In particular, u — (u\ 1 1*2,113) has weights (5,3,1) in the (2, 7)-case and 
weights (5, 2, 1) in the (3, 4)-case. 

We now consider functions that are periodic with respect to the lattice A. 

Definition 2.2. Let 9Jt(u) be a meromorphic function ofu e C 9 . Then 9Jt is a standard Abelian function 
associated with C if it has poles only along k _1 (0[ 9_1 ]), and satisfies 

M(u +£) = Tt(u) for all t E A. (19) 

Note the comparison with equation ([T]) and that the period matrices play the role of the scalar periods in the 
elliptic case. We will define generalisations of the Weierstrass p- function which satisfy equation (fl~9| , defined 
using the quasi-periodic function defined below. 

First, let 8 = S'oj' + 8"oj" be the Riemann constant with base point 00. Then g„ is the theta charac- 
teristic presenting the Riemann constant for the curve C with respect to the base point 00 and generators 
{ay, f3j} of H\(C,Z). For any ((2, 7)-curve, we have a classical choice of {ctj,(3j} and (see [3T] or [TT]) we 
have 

<5' = [|l|r, ^[Hf (20) 
For any (3,4)-curve, Shiga computed in [25| that with his choice of {ex,-, /?-,-} we have 

8' = [0 I 0] T , 8" = [0 \ 0] T . (21) 

Definition 2.3. The Kleinian cr-function associated to a genus three (n, s)-curve may be defined using a 
multivariate 9 -function with characteristic 8 as 

o-(u) = cc^p(^ut 1 '{lj , )~ 1 u t ) ■6[8}((lj')' 1 u t \ (aZ^V)- 
= cexp (iit?7'(a/) -1 it T ) 



X 

TTlE 



J2 exp 



2m{\{m + d'Hw')" V'(m + $') + ( m + 8') T {{lo 1 )- 1 ^ + 8")} 



The constant c is dependent upon the curve parameters and the basis of cycles and is fixed later, following 
Lemma\T 
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We now summarise the key properties of the cr-function. Sec [S] or [22] for the construction of the 
cr-function to satisfy these properties. 

For any point u G C 3 we denote by v! and u" the vectors in R 3 such that 

u = u'u' + u"u)". (22) 

Therefore a point £ G A is written as 

£ = I! J + £"lu" g A, £', £" e Z 3 . (23) 
For u, v G C 3 and £ € A, define u) and x(^) as follows: 

w) = u(t?'r?' + v"rj") T , 
Xifi = exp [2ri{l'{5") T - r'(<5') T + Wf}] • 
Lemma 2.4. TTie a -function has the following properties. 

• It is an entire function on C 3 . 

• It has zeros of order one along the set k~ 1 (Q^). Further, cr(u) ^ outside this set. 

• For all u G C 3 ,£ G A the function a(u) has the following quasi-periodicity property: 

£ 



a(u + £)= x(0 cxp 



Liu, 



(24) 



• It /las definite parity given by a([—l]u) = ±a(u) in the (2,7) and (3,4)-cases respectively. 

Proof. The function is clearly entire from the definition, while the quasi-periodicity and zeros of the function 
are classical results, (see [3]), that are fundamental to the definition of the function. They both follow from 
the properties of the multivariate ^-function. The parity property follows from Proposition 4(iv) in |22j which 
states that for a general (n, s)-curve, 

a(-u) = (-l)^™ 2 - 1 ^ 2 - 1 ^). (25) 

□ 

We next define p- functions using an analogy of equation ([5]). Since there is more than one variable we 
need to be clear which we differentiate with respect to. We actually define several multivariate p-functions 
and introduce the following notation. 

Definition 2.5. Define m-index Kleinian p-functions, (where m>2) as 

i \ d d d r . x1 r , 

Pii,i2,...,i m { u ) = ~q^~q^~ ■ ■ -~du~ g L a WJ' *!<•••< hn e{l,...,g\. 

The m-index p-functions are meromorphic with poles of order m when a(u) = 0. We can check that they 
satisfy equation (|19[) and hence they are Abelian. The m-index p-functions have definite parity with respect 
to the change of variables u — > [— l]u and are odd if m is odd and even if m is even. Note that the ordering 
of the indices is irrelevant and so for simplicity we always order in ascending values. 

When the (n, s)-curve is chosen to be the classic elliptic curve then the Kleinian cr-function coincides with 
the classic cr-function and the sole 2-index p-function coincides with the Weierstrass p-function. The only 
difference is the notation with 

Pii(u) = p(u), pm(w) = p'(u), pim(u) = p"(u). 

Clearly, as the genus of the curve increases so do the number of associated p-functions. For the curves we 
study we set g = 3 to leave six 2-index p-functions, ten 3-index p-functions etc. By considering Definition [53] 
we see that the weight of the p-functions is the negative of the sum of the weights of the variables indicated 
by the indices. So although the two curves have the same number of associated p-functions, they are of 
different weights arising from the different basis of holomorphic differentials. We will find in the following 
sections that this leads to variations in the theory of the two sets of functions. 

We now introduce a final result detailing how the functions can be expressed using series expansions. 
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Lemma 2.6. The Taylor series expansion ofa(u) about the origin may be written as follows: For each type 
of cyclic genus three curve, there are constants K 2t 7 and K 3A depending only on the Xj and {oj, flj} such 
that 



a(u) = < 



K 2 ,7 ■ (sWzjiu) + ^2 ^e+2k(u, A)) for a cyclic (2, 7)-curve 1 

k ^ (26) 
K 3A ■ (^SW 3A {u) + ^ C 5+3k (u, X)j for a cyclic (3, A)-curve. 



k=l 



Here SW n _ s is the Schur- Weierstrass polynomial generated by (n, s) and each C m is a polynomial composed 
of products of monomials in u of weight m multiplied by monomials in the curve parameters of weight (6 — m) 
and (5 — m) in the (2,7) and (3,4)-cases respectively. 

A similar result may be stated for the non-cyclic curves using fij instead of , but the calculations involved 
in deriving such an expansion are computationally much more difficult. 

Proof. We refer the reader to [32] for a proof of the relationship between the cr-function and the Schur- 
Weierstrass polynomials and note that this was first discussed in [5]. We see that the remainder of the 
expansion must depend on the curve parameters and split it up into the using the weight properties. 
Each Cfe must be finite since the number of possible terms with the prescribed weight properties is finite. 

□ 

Expansions of this type were first introduced in [7] in relation to the study of the Benney equations. 
Since then they have been an integral tool in the investigation of Abelian functions. Recently, computational 
techniques based on the weight properties have been used to derive much larger expansions and we refer the 
reader to [T^ and [TB] for details of how to construct and use these expansions. 

We fix c in Definition 12.31 to be the value that makes the constant one in the above lemma. Some other 
authors working in this area may use a different constant and in general these choices are not equivalent. 
However, the constant can be seen to cancel in the definition of the p-functions, leaving most results indepen- 
dent of c. Note that this choice of c ensures that the Kleinian er-function matches the Weierstrass cr-function 
when the (n, s)-curve is chosen to be the classic elliptic curve. 

The connection with the Schur- Weierstrass polynomials also allows us to determine the weight of the 
cr-function as (l/24)(n 2 — l)(s 2 — 1). In the (2, 7)-case this gives a(u) weight 6 while in the (3,4)-case it has 
weight 5. The respective Schur- Weierstrass polynomials are, 

SW 2 , 7 = ±u 6 3 - |u|u 2 - u\ + U 3 ttl, (27) 

SW 3A = sjug-usua + ui. (28) 

The a- function expansion has been calculated up to and including the polynomial C 3 % in the (2, 7)-case and 
C35 in the (3, 4)-case. These expansions are available from the authors and should also be available from the 
journal in the Supplementary Materials to the paper. 



3 Bases for the vector spaces of Abelian functions 

3.1 General theory 

We classify the Abelian functions according to their pole structure. We denote by r(j, 0(mG' fc ')) the vector 
space of Abelian functions defined upon J which have poles of order at most m, occurring only on the fcth 
strata, 0^. The case where k = g — 1 is of interest since all the Abelian functions we deal with have poles 
occurring here, on the 0-divisor. A key problem in the theory of Abelian functions is the generation of bases 
for these vector spaces. 

Note that the dimension of the space T(J, 0(m6' r1 ')) is m 9 by the Riemann-Roch theorem for Abelian 
varieties. (See for example [20] page 99.) Recall that the m-index p-functions all have poles of order m. We 
see that the number of m-index p-functions associated with a genus g curve is 

(g + m-l)\/[ml(g-l)\], 

which will not grow as fast as the dimension of the space. We hence need to identify a wider class of Abelian 
functions than the p-functions in order to construct such bases. 
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Recall that an entire Abelian function must be a constant and that there is no Abclian function with a 
single pole of order one. Hence those Abelian functions with poles of order two are the simplest and so are 
often referred to as fundamental Abelian functions. The basis problem has been solved in general for such 
functions, through the inclusion of the following extra class of Abclian functions. 

Definition 3.1. Define the operator T>i as below. This is now known as Hirota's bilinear operator, although 
it was used much earlier by Baker in J2J/. 

d d 
T>i = . 

dui dvi 

Then an alternative, equivalent definition of the 2-index p-functions is given by 

1 



pijiu) 



T ViT>j<j{u)a{v) 



i<je{l,...,g}. 



2a(u) 2 ' 

We extend this approach to define the m-index Q-functions, for m even, by 

(-1) 



Qii ■ 



where i% < ■ ■ • < i m £ {1, . . . , g}. 



2a(u) 



.T> im a(u)a(v) 



(29) 



(30) 



The m-index Q-functions are Abelian functions with poles of order two when o~(u) = 0. Note that if you 
were to apply the definition with m odd then the functions would be identically zero. A 4-index Q-function 
was originally used by Baker with the generic 4-index Q-functions introduced when research first started 
on the trigonal curves. (In the literature they are just defined as Q-functions). The definition for m-index 
Q-functions above was developed in jTS] as increasing classes are required to deal with cases of higher genus, 
(see also [16]). In this paper we only need to use 4-index Q-functions, which in [13] were shown to satisfy, 



Qijki — Pijki — 2pijpke — 2pikpje — 2pupjk- 



(31) 



Similar expressions have been found for the higher index Q-functions, (see |18j). 

When considering the vector spaces for Abelian functions with poles of higher order, a natural place to 
look for extra functions is in the derivatives of the functions with lower order poles. Note that while the 
derivatives of m-index p-functions are (m + l)-index p-functions, the same is not true for the Q-functions. 
For brevity we adopt the notation 

d 

diQ l3 ki{u) = -g^-Qijki(u) 

and similarly for other functions. As discussed in [53], the derivatives of existing basis functions will not 
be sufficient to find successive bases in the case where the theta divisor has singular points. In [T3] the 
authors introduced the following new class of functions to overcome this problem for the three pole basis in 
the (3,4)-case. 

Definition 3.2. Consider the matrix, 



[Pij] 



3x3 



pll pl2 pl3 
P21 P22 P23 
P31 P32 P33 



Pll Pl2 Pl3 
Pl2 P22 P23 
Pl3 P23 P33 



Then the function is defined to be the minor of [pij]3xi 



For example, 



,[12] 



Pl2 
Pl3 



P23 
P33 



P12P33 - P23pl3- 



So each of these functions will be a sum of two pairs of products of 2-index p-functions. Although each 
product will have poles of order four, we can check by substituting with Definition 12.51 that these cancel so 
that the function has poles of order three. Although this class solved the three pole basis problem in the 
(3, 4)-case, it was not sufficient to complete the corresponding basis in the (2, 7)-case. We explicitly address 
this and similar problems for the four pole vector spaces below. A general method to generate new Abelian 
functions is currently being developed with some recent new results in [TT] . 
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3.2 The general (2, 7)-curve 

In this section we treat the most general (2, 7)-curve, namely the curve defined by equation @. A well known 
basis for T(j, O(20^)) in the (2, 7)-case is 

CI © Cpn © Cpi2 © Cpi 3 © Cp 22 © Cp 23 © Cp 33 © CA, (32) 

where 

A = pnp33 - Pl2p23 ~ Pl3 + Pl3p22- (33) 

The function A has poles of order three in general, (we can this check using Definition 12.51) . However, these 
cancel to order two in the (2, 7)-case only, (which we can check using the cr-expansion). The A- function was 
introduced by Baker in [4] although he did not use it explicitly to form such bases. More recently the A- 
function has been used in [TJ in a covariant analogue to the theory of hyperelliptic p-functions. It is possible 
to rewrite the theory to match the general approach suggested above, by replacing A with a Q-function of 
weight —12. However, the use of A is advantageous since it allows the theory to be completely realised in 
terms of 2 and 3-index p-functions, with all higher index p-functions given recursively in these. 

The construction of r( J, 0(38 1 2 ')) has recently been studied by Nakayashiki in [23]. The author was able 
to enumerate in detail all the terms with the exception of the final one, which was labeled F 3 , (see [33] page 
27 onwards). It was shown that 

1 + (lower degree terms) 



and that it satisfies certain differential equations involving power series of the u t . In Theorem 13.31 below we 
give a new explicit form for this term, which we rename as T, in terms of p-functions. 

Theorem 3.3. The basis for T(j, O{30 [2] )) is given by 

d © Cp m © Cpua © Cpn 3 © Cpi 22 © Cp 123 © Cp 133 

© Cp 222 © Cp 223 © Cp 23 3 © CP333 © C<9iA © C<9 2 A (34) 
© C<9 3 A © CpM © Cpl 12 l © Cpl 13 l © Cpl 23 l © Cpl 33 ! © CT, 

where 

T = P222 - 4 P22 - Q2222P22 = 2p2 2 + p 222 - p 22 p 2222 . (35) 

Proof. The dimension of the space is 3 9 = 3 3 = 27 by the Riemann-Roch theorem for Abelian varieties. 
All the selected elements belong to the space and we can easily check their linear independence explicitly in 
Maple using the cr-expansion. 

To actually construct the basis we started by including the 8 functions from basis (|32|) for the functions 
with poles of order at most two. We then know that the remaining entries must have poles of order three. 
We start by looking for entries from the set of derivatives of the basis ([32]) . We test at decreasing weight 
levels and look to see whether these functions can be written as a linear combination upon substitution of 
the series expansions. (A more detailed discussion of the algorithm used for this is available in [T5].) 

Note that while this theorem holds for the general (2, 7)-curve in equation ([9]), we need only use the series 
expansions associated to the cyclic (2,7)-curve in equation pip . This is because if an element cannot be 
expressed using the basis with the restriction on the parameters, then neither will it be expressible with the 
wider set of parameters. Further, we only need to use sufficient expansion to give non-zero evaluations of the 
functions considered in order to check whether they are linearly independent. 

After examining all these functions we find that 21 basis elements have been identified. We identify a 
further 5 by considering the class given in Definition 13.21 Note that unlike the (3,4)-case, here 

P [131 = P12P23 ~ P22P13 and p [22] = P11P33 - PI3 

have the same weight, (—12). Further, they are linearly dependent and so only one can be included in the 
basis, (either may be chosen). 
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To find the final function the following new class of functions is considered. 

Tijklmn = PijkPlmn ~ ^PijPklPmn ~ \pijPkmPln ~ \pijPknPlm ~ ^PikPjlPmn (36) 

2 2 2 1 1 

— ^PikPjmPln ~ 3 PikPjnPlm ~ 3 PilPjkPmn + 3 PilPjmPkn + 3 PilPjnPkm 

2 112 11 

— zPimPjkPln + 3 PimPjlPkn + 3 PimPjnPkl ~ 3 PinPjkPlm + 3 PinPjlPkm + 3 PinPjmPkl 

gQijklPmn 3 QijkmPln 3 QijknPlm © 3 QijlmPkn © 3 QijlnPkm © 3 QijmnPkl 
© 3 QiklmPjn © 3 QiklnPjm © 3Q ikmnPjl 3 QilmnPjk © 3 Q jklmpin © "^Q jklnpim 
© 3 QjkmnPil 3 QjlmnPik 3 QklmnPij- 

Substituting with Definition 12.51 shows these functions to have poles of order at most three. We stress that 
this is the case for the T-functions associated to any curve, and not just the genus three case we consider 
here. The functions were derived through an attempt to match in general, the poles of a quadratic term in the 
3-index p-functions with a polynomial in 2-index p-functions. Although this is not possible, we can match 
by including Q-functions as well. A general polynomial of this type was constructed and the coefficients then 
determined as above to ensure the poles of order greater than three vanish. A method to generate new basis 
functions is currently being developed with some recent results published in [T7] . 

The new class of functions is examined at decreasing weight levels. We find that all can expressed as a 
linear combination of existing basis entries except for those at weight —18. We can choose any T-function 
at weight —18 to complete the basis and use the simplest of these, 7222222 = T as given in equation (|35]l. 

□ 

Theorem 3.4. The basis for T( J, 0(4G>[ 2 ])) is given by 

© © Cpiui © Cpma © Cp m3 

© Cpii 33 © Cpi222 © Cpi223 
© Cp2222 © Cp2223 © Cp 2 233 

© O13A © O23A © C<9 33 A 

© C<9 U A © Oipl"! © 2 p [11] 

© C9 3 p [12] © Cdipl 13 ! © 2 pl 13 l 

© Cd 2 p [23] © C9ipl 33 l © 2 pl 33 l 

© Cd 3 T © CG, 

where G = P22222222 - 140p§ 222 . 

Proof. We follow the proof of Theorem 13.31 This time the dimension is 4 9 = 4 3 = 64 and we find that we 
can identify 63 functions using the basis (|34j) and its derivatives. 

To find the final basis function the following new class of Abclian functions was examined. 

Qijklmnop — Pijklmnop ^(pijnopklmp © PijnpPklmo © PijlpPkmno © PijmnPklop © Pijmopklnp 
© PijmpPklno © PijkmPlnop © PijknPlmop © PijkoPlmnp © PijkpPlmno © PijlmPknop 
© PijlnPkmop © PijloPkmnp © PijklPmnop © PijopPklmn + P iklm^Pjnop H~ ftikln'Pjmop 
fpiklofPjmnp $iklp$jmno + p t kmnfPjlop $ikmo < (Pjlrip 'Pikmpfpjlno fpiknofpjhnp 
fPiknpfPjlmo fPikopfPjlmn + p ilmnPjkop © PilmoPjknp © PilmpPjkno © PilnoPjkmp 
© PilnpPjkmo © PilopPjkmn © PimnoPjklp © PimnpPjklo © PimopPjkln © PinopPjklm) • 

The C?-functions associated with any curve have poles of order at most four (using Definition 12. 5j) and were 
derived by matching the higher order poles in an 8-index p-function using arbitrary sum of quadratic terms 
in the 4-index p-functions. Examining at decreasing weight levels we see that we need to include a function 
at weight —24 and choose £22222222 = G as given in the theorem to complete the basis. 

□ 

We note that these bases are not unique. As discussed at the start of the section, we could replace the 
A-function with a Q-function if desired. Similarly, in the theorems above we could have used any T-function 
at weight —18 and any (/-function at weight —24 to complete the 3 and 4-pole bases respectively. However, 
we can conclude that the weight structure of each basis is unique. 



© Cpn 2 2 

© Cpi233 

© Cp2333 

© C5i 2 A 

© ca 3 p [11] 

© C<9 3 p [13] 
© CdiT 



© Cpn23 

© CP1333 

© CP3333 

© ca 22 A 

© Cdipl 12 ! 

© Cdipl 23 ] 

© Cd 2 T 
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Lemma 3.5. Consider a basis for T(j,0(kQ^ 9 associated to any (n, s)- curve. Scale the functions so 

that they do not depend on any curve constants. Then the set of weights of the functions is unique. 

For example, the basis (|32p has weights [— 12, —10, —8, —6, —6, —4, —2, 0]. While the functions for this 
basis may be changed, this set of weights can not. 

Proof. Suppose that we could replace a basis entry with another of different weight. The replacement function 
would have to be linearly independent of the other basis entries, and hence a constant multiple of the original 
function. This would lead to equations that are inhomogcneous in the Sato weights. 

□ 

We can also determine a minimal bound of these weights, as detailed in the following lemma. 

Lemma 3.6. The basis for the space T(J, 0(fc0[ ff-1 l)) associated to an (ri,s)-curve will contain functions 
with weights no lower than — fcwt(cr). 

Proof. Such basis functions may be written as a rational function in a(u) with overall denominator a(u) k . 
The simplest possible numerator is a constant, which gives a function with weight — fcwt(er). Any other 
numerator would depend on u and so have a higher weight. 

□ 

Lemma 13.61 stops us from testing those functions at a lower weight thus drastically reducing the amount 
of computation required to find such bases. These lower weight functions must be expressible as a linear 
combination of the basis functions in which every term depends on the curve parameters. 

We conjecture that every such basis may be evaluated using polynomials of p-functions. This is known 
to be true for the curves considered here (see [5] and [TU]) and should follow similarly for all (n, s)-curves. 
All such polynomials have negative weight except the constant function which will take the maximal weight 
of zero. Together with Lemma 13.61 this would restrict basis entries to the weight range — fcwt(cr), ... ,0. 

We note that in the (2, 7)-case the weight of cr{u) is +6 and the functions A, T and G have weight 12, 18 
and 24 respectively. They are hence the minimal weight entries in their respective bases. By considering 
the rational limit of these functions, where o~(u) = SW2^(u), we can check that these functions each have 
leading term of a constant over a(u) k with k = 2, 3, 4 respectively. 

It is clear that functions of the form (const + 0(ui))/a k play a special role in completing the construction 
of r( J, 0(fc9[ 9_1 l )) . Such functions are straightforward to construct using the rational limit and the method 
of undetermined coefficients. For example, the following function evaluates to a constant over er 5 and we 
postulate that such a function will be essential for the 5-pole basis. 

12p22 - 8pi 2 P2222 + 6pp22p222 + P22P2222 ~ P2222P222- 

3.3 The general (3,4)-curve 

Here we treat the most general (3,4)-curve, namely the curve defined by equation ([T0| . Lemma 8.1 in [13] 
identified the basis for T( J, C(29 [21 )) in the (3,4)-case as 

CI © Cp n © Cpi2 © Cpi 3 © Cp 22 © Cp 23 © Cp 33 © CQ1333, (37) 
and the basis for T(J, 0(39 M)) as 

(EZD © Cpm © Cpua © Cpii 3 © Cpi 22 

© Cpi 33 © Cp222 © Cp 2 23 © Cp 23 3 
© CdiQi333 © C9 2 Ql333 © Cd 3 Ql333 © Cp^ 

© Cpl 13 l © Cpl 22 l © Cpl 23 ! © Cpl 33 l 

Note that here the functions pl 13 l and p' 22 ' have weight —10 and —12 respectively. They were hence inde- 
pendent meaning this class was sufficient to complete the basis and no 7~-fuiictions were required. However, 
we find that upon proceeding to the four pole basis we will have to define an extra class. 



-'123 



^P333 

Cpl 12 ! 



(38) 



10 



Theorem 3.7. The basis for T( J, C(46 [21 )) is given by 

(USD e Cpim © Cp m2 © Cp m3 

© CpH23 © CpH33 © Cpl222 

© Cpi233 © Cpi 333 © Cp2222 

© Cp2233 © Cp2333 © Cp 33 33 

© C<9i<9 2 Qi333 © C9ia 3 Ql333 © C5 2 9 2 Ql 3 33 

© Cd 3 d 3 Q 1333 © Cd lP W © Oipl 12 ! 

© cd lP w © o lP [ 23 i © o lP i 33 ] 
© ca 2 p [12] © o 2 p [13] © o 2 p[ 22 i 
© o 2 p [33] © cd 3 p [11] © C9 3 p [121 

© CF, 

where 

F = Pllp22p33 - P11P23 ~ P12P33 + 2pl2pl3p23 ~ P? 3 P22- (40) 

Proof. We follow the proof of Theorem 13.31 This time the dimension is 4 9 = 4 3 = 64 and we find that we 
can identify 63 functions using the basis (|38[) and its derivatives. 

To find the final basis function the following new class of Abclian functions was examined. 

Fijklmn = PijPklPmn ~ PijPknPlm ~ Pilpjkpmn + PilPjnPkm 

~t~ PimPjlPkn PimPjnPkl ~t~ PinPjkPlm PinPjlPkm- 

The J-"-functions associated with any curve can be seen to have poles of order at most four by substituting 
with Definition 12.51 and were derived by matching the higher order poles in an arbitrary sum of cubic terms 
in the 2-index p-functions. Examining at decreasing weight levels we see that we need to include a function 
at weight —16 and choose .F112233 = F as given in equation (|4TJ|) to complete the basis. 

□ 

This new basis allows us to derive the new addition formula in Theorem 15.31 

The problem of determining such bases is of importance and is a barrier to further development of the 
theory. The construction of the T, jF and C?-functions in this section represent great progress towards a 
general solution of the problem. We have used these functions, along with another new class to solve the 2 
and 3-polc basis problems for the (2,9)-curve (genus four) and present these in Appendix [Al 

4 Differential Equations 

The p-functions associated to (n, s)-curves satisfy a variety of differential equations which we review in this 
section. Some of these differential equations can be found occurring naturally in areas of mathematical 
physics, so the p-functions can be used to give solutions to a variety of important problems. In this section 
we consider the three main classes of differential equations and compare the explicit equations for the two 
genus three cases. Note that although the functions associated to the (2, 7) and (3, 4)-curves are notationally 
the same, they behave differently and satisfy different equations. This is most apparent from the different 
weights assigned to the functions, summarised for the fundamental functions in the table below. 





P11 


P12 


P13 


P22 


P23 


P33 


A 


Q1333 


(2, 7)-case 


-10 


-8 


-6 


-6 


-4 


-2 


-12 




(3, 4)-case 


-10 


-7 


-6 


-4 


-3 


-2 




-8 



The sets of differential equations in this paper are presented in decreasing weight order as indicated by 
the bold number in brackets. They have all been made available in text files in the supplementary material. 
Note that they refer to the functions associated with cyclic curves, but that relations for the general curves 
can be derived in a similar fashion at a greater computational cost. 

We first consider the set of differential equations to express the 4-index p-functions. These J^-index rela- 
tions are the generalisation of equation ((4]) from the elliptic case. We aim to express each 4-index p-function 
as a degree two polynomial in the 2-index p-functions, in comparison with equation ((4]). Then, through 
differentiation and manipulation of this set, we could express all higher index p-functions as polynomials in 



CpH22 

Cpl223 

Cp2223 
C<9i<9i(2i333 
C<92<9 3 Qi333 

Cd lP W (39) 

o 2 pi n ] 
ca 2 p[ 23 i 

C9 3 P [221 
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the 2 and 3-index p-functions, analogous to the elliptic case. A complete set of such relations can be obtained 
for the (2, 7)-case as first presented by Buchstaber, Enolskii and Leykin in [5]. 

(-4) p3333 = 4p23 + 4p 33 A 6 + 2A 5 + 6P33 (41) 
( — 6) p2333 = 6pi 3 - 2p 22 + 4p 23 A 6 + 6p23p33 (42) 

( — 18) piu2 = 6pnpi2 - 4A 5 A - 2p22Ai + 6pi 3 Ai - 8p 2 3A + 4pi 2 A 2 
(-20) pun = 6p 41 + 2A 3 Ai - 8A 4 A + 4p u A 2 + 4p i2 Ai - 12p 22 A + 16pi 3 A 

The complete set is displayed in Appendix [B] There are various ways to derive such relations, with a 
recent survey given in [15]. A constructive way is to consider the basis for r( J, 0(20 1 2 ')) and the set of 
4-indcx Q-functions. Each Q-function has poles of order at most two and so belongs to the vector space 
r(j,e>(2el 2 ])). Hence each Q-function can be expressed as a linear combination of basis entries. (The 
explicit linear combination can be identified using the tr-expansion as discussed in |18|). In the (2, 7)-case we 
have basis (|32|) containing 2-indcx p-functions and A. Hence we can substitute for the Q-functions and A 
using equations (|3"Tj) and (j3"3"j) to leave the desired set of 4-index relations. 

The corresponding set of equations for the (3,4)-case was derived in Lemma 6.1 of [13] and are also 
presented for comparison in Appendix [Bl 

(-4) p 3333 = -3p 22 + 6p 33 
(—5) P2333 = 3p 33 A 3 + 6p 23 p 33 

( — 17) pni2 = -Qi 333 Ai + 6p 33 A 3 A + 6pnpi2 

(—20) pun = -4<2i 3 33Ao - 3p 33 Ai 2 -I- 12p 33 A 2 A + 6pli 

Note that here the function Qi 333 is used in some of the expressions. The reason for its inclusion can be 
explained by considering the constructive method discussed above. Since Qi 333 was used in basis (|37[) for 
r(j,e>(2ei 2 i)) it win appear in the expressions following from this basis. The ability to construct 4-indcx 
relations using only 2-index p-functions as in the (2, 7)-case is a feature that appears to be unique to the 
hyperelliptic cases. We have explicitly checked that such relations cannot be achieved at certain weights in 
the (3,4)-case. A more appropriate definition for 1^-index relations seems to be a set that expresses all 
the 4-index p-functions using a degree two polynomial in the fundamental basis functions. (Note that the 
Q-functions used will only need to appear linearly.) 

It is natural to next consider a set of differential equations to generalise equation ([3]). Such a set should 
give expressions for the product of two 3-index p-functions and so we refer to them as quadratic 3-index 
relations. The natural generalisation would express each product as a degree three polynomial in 2-index 
p-functions, but as in the previous case, this appears to only be possible for the hyperelliptic functions. We 
have again explicitly checked that such relations do not exist at certain weights in the (3, 4)-case and so 
propose the modified definition of quadratic 3-index relations to be a set of differential equations that 
expresses all the products of 3-indcx p-functions using a degree three polynomial in the fundamental basis 
functions. 

The quadratic 3-index relations for the (2, 7)-case were considered in [8], but a complete set was not 
presented. In [T] the corresponding relations for functions associated to covariant curves have been considered, 
but again a complete set is not directly obtainable from the results published there. 

Theorem 4.1. The quadratic 3-index relations associated to the cyclic (2,7)-curve start with those below, 
with the complete set as displayed in AvvendixXCk 



(- 


-6) 


2 

P333 


= 4p| 3 + 4A 4 + 4A 5 p 33 + 4A 6 p 2 3 - 4p i3 + 4p 22 + 4p 33 p 23 (43) 


(- 


-8) 


P233P333 


= 2(2p 23 p 33 + A 3 + A 5 p 23 + 2A 6 p 33 p 23 + P12 + 2p 33 pi 3 - p 33 p 2 2 + p| 3 ) 


(- 


-10) 


2 

P233 


= 8p 23 pi 3 - 4p 23 p 22 + 4pn + 4A 2 + 4p|jp 33 + 4A 6 pl 3 


(- 


-10) 


P223P333 


= ~4p 33 pi2 - 2p 23 pi 3 + 4p 23 p22 + 4A 6 p 33 pi 3 + 2A 5 p 33 p 23 + 4A 4 p 23 








- 2pn + 2A 5 pi 3 - 2A 3 p 33 + 2p 22 p| 3 + 2p| 3 p 33 


(- 


-10) 


P133P333 


= -2p 33 pi2 + 2p 23 pi 3 + 4A 6 p 33 pi 3 - 2pn + 2A 5 pi 3 + 4pi 3 p^ 3 
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Proof. Once again, these relations can be derived through a variety of methods as discussed in |12j . They can 
again be found constructively using the c-expansion, although in this case there is no simple linear algebra 
result to dictate that such relations exist. Instead we may just search for them using arbitrary polynomials 
of 2-index p-functions. The computations involved can be heavy and so simplifications are made by ensuring 
the polynomials are homogeneous in weight. Additionally, it is often possible to allocate the cubic terms to 
cancel higher order poles using Definition 12.51 at a low computational cost. 

□ 



Lemma 4.2. Theorem\^.l\ may be written concisely as a determinantal formula, 



(l T Ak)(l' T Ak') 



H 


V 


k' 


l T 








k T 









(44) 



using the matrices 



A = 





P333 
— P233 
P223 — Pl33 
"P222 + 2pi23 



Pl33 
-P123 
Pl22 — Pll3 



H = 



4A 
2Ai 
-2pn 

-2pl2 

-2pi3 



2Ai 
4A 2 + 4 
2A 3 

"2p22 +4 
-2p23 



4P11 

2pi2 



"P333 P233 
-Pl33 


PU3 
"Pll2 

-2pn 
2A 3 + 2pi2 
4A 4 + 4p 22 - 4pi 3 
2A 5 + 2p23 

-2p33 



;i3 



P223 + Pl33 
Pl23 
-Pll3 



pill 

-2pl2 
-2p22 + 4pi 3 

2A 5 + 2p 23 
4A 6 + 4p 33 
2 



P222 ~ 2pl23 
-Pl22 + Pll3 
P112 
"Pill 





-2pl3 
"2p23 
"2p33 

2 




and arbitrary vectors l,k,l',k' of dimension 5. 
using appropriate choices of the vectors. 



Each of the 55 relations in Appendix [0 may be obtained 



Proof. The formula (|44[) may be verified directly by expanding and substituting with the relations of Theorem 
14.11 To derive the individual relations from the determinantal formula one must make suitable choices of the 
vectors Z, k, I' , k' with only one entry in each vector non-zero. For example, setting 

k = k T = (1,0,0,0,0), l = l T = (0,1,0,0,0) 

into equation (|44l) leaves relation f|43[) for p| 33 . One restriction on this approach is that some equations 
must be derived before others. For example, the equation for pf 33 must be derived before the equation for 
P133P223 as the choice of vectors to give P133P223 on the left hand side will introduce pi 33 as well. However, 
when this is the case it is always possible to derive one of the other quadratic term independently and then 
substitute for it to find the other. Hence, all the 55 individual relations may be derived from the formula. 

□ 

Lemma 14.21 was inspired by the results of Athorne in [1] on covariant hyperelliptic curves and the corre- 
sponding p-functions. These curves and functions belong to generic families permuted under an SI2 action 
which can be easily mapped to the standard curves and functions considered here. Recently, the covariant 
result corresponding to Lemma 14.21 was developed in [2J. 

The corresponding quadratic relations for the (3, 4)-case were first considered in [13] . but again a complete 
set has not been available until now. These relations include the basis function Q1333 which occurs only 
linearly, or multiplied by a single 2-index p-function. A determinantal version of these equations, one similar 
to equation (|44|) . has not been identified. One may follow from development of the corresponding covariant 
theory to trigonal curves. 
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Theorem 4.3. The quadratic 3-index relations associated to the cyclic (3,A)-curve starts with those below 
with the complete set as given in AvvendixWi 



-6) 


2 

P 333 


= 4 P33 + 4 Pl3 + P23 _ 4p22p33 




-7) 


P233P333 


= 4p23pl 3 - 2 Pi2 - P22P23 + 2A 3 p§ 3 




-8) 


P233 


= 4p| 3 p33 + P22 + 4A 3 p 23 p33 + 4p 33 A 2 - §Ql 333 




-8) 


P223P333 


= 2 P23P33 + 2p22p§ 3 ~ 2p| 2 + A 3 p 23 p 33 + 4p i3 p 33 + |C 


?1333 


-9) 


P223p233 


= 2pi 3 A 3 + 2Al + 4pi 3 p23 + 2p 23 p22p33 + 2p 2 ! 3 

+ P22P33A3 + 2p^ 3 A 3 + 2p 23 A 2 




-9) 


P222P333 


= 6p 23 p22p33 - 4pi2p33 ~ 2p^ 3 + 4p 22 p 33 A 3 - 8pi 3 p 23 


- P23 A 3 



'13^3 



Proof. These relations can be derived using a variety of methods as in Theorem 14.11 



□ 



The final set of differential equations that we consider here are a set bilinear in the 2 and 3-index p- 
functions. Due to the parity properties of the p-functions we know that these bilinear relations cannot 
contain any constant terms, or terms dependent only on the 2-index p-functions. There is no analogue of 
these relations in the genus one case. Three of the (2, 7) bilinear relations displayed below were derived in 
[5] as part of a larger computation, but the complete set below is a new result. 

Theorem 4.4. Every bilinear relation associated with the cyclic (2,7)-curve may be given as a linear com- 
bination of the 24 below. 



(-7 
(-9 
(-9 
-11 
-11 

-11 
-13 

-13 
-13 
-13 
-15 

-15 
-15 

-15 
-17 

-17 

-17 
-17 



= P233P33 + P223 ~ P333P23 - P133 

= P133P33 + Pl23 - P333P13 

= -2p 33 p 223 + P22P333 + P23P233 + A5P233 ~ 4p i23 + 2A 6 pi 33 + p 2 22 ~ 2A 6 p22 3 
= -pl 33 p23 - Pll3 + P233P13 

= p 33 p222 — 2p22p233 + P23p223 + A.5p223 - Pl22 + 2A 6 pi 23 + 3pn 3 — 2Aspl 33 

+ P333A3 - 2p 233 A4 

= — P133P23 - 2p 33 p i23 - 3pH3 + P122 - 2A 6 pi 23 + pl2p 333 + 2p 233 pi 3 + A5P133 

= P112 - P222P23 + P22P223 + P233A3 + 2pi 33 pi 3 - 2p 22 pi 33 + 2pi 2 p 233 

- 2pnp 333 - 2p 333 A 2 - 2A 6 pn 3 

= -P123P23 ~ P133P13 + P13P223 + P33P113 

= -p22pl33 - 2pi23p23 + Pll2 — 2A 6 pn 3 + pl2p233 + 2pi3p223 

= P122P33 - 2p22pl33 + A 5 pl23 + P133P13 + P13P223 - P11P333 - 2P133A4 

= — A5P113 + 2pi 33 A 3 + 2pi 2 p 223 + P22P123 — 3p 333 Ai + 2pi 23 pi 3 

- 2p222pl3 - P122P23 - P112P33 - P11P233 

= — A5pn 3 + pi 33 A 3 + P12P223 — P333A1 + pm — P122P23 

= — P133A3 + 4pi 23 pi 3 - 2pi 2 pi 33 — pllp233 - P23P113 + Pl22p23 

+ Pl2p223 - p222pl3 ~ p22pl23 
= P133A3 + P12P133 + P11P233 - P23P113 - P112P33 - A5P113 
= -A5P112 + P12P222 - P22P122 - P233A1 - 4pi 3 pn 3 + 2pi 3 p i22 + 2p 22 pn 3 

- 2p 23 pn 2 - 2p i2 pi 23 + 2p 33 pm + 2pnpi 33 + 2pi 33 A 2 + 2A 6 pm 

= -P22P113 + P11P223 + 3pi2pi23 - P233A1 + 2pi 33 A 2 - 4p 333 A + pi 23 A 3 

+ 2pi 3 pn 3 - 2pi 3 p i22 - p2 3 pii2 - P33P111 - P11P133 - 2pn 3 A 4 
= P22P113 + 2pi2pi23 - P233A1 + 2pi 33 A 2 - 4p 333 A - 2pi 3 p i22 - p2 3 pii2 
= -P233A1 - 3pi 3 pn 3 + 3pnpi33 + 2pi 33 A 2 + 2p 333 A - P23P112 - P12P123 
+ P13P122 + P22P113 
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( — 19) = -2p22pll2 + pllp222 + P12P122 - Aip223 + 2A 2 pi23 - 4p233Ao + P122A3 

+ A5P111 - 2A4P112 - 2pii 3 A 3 + 3pi 33 Ai 
(—19) = P133A1 + P12P113 - 2P233A0 - P13P112 

( — 19) = -3pi 33 Ai + A1P223 - P12P113 - 2pnpi23 - 2A 2 pi23 + 4p233Ao + P113A3 

+ 2pi3pH2 + p23plll 

(—21) = -4Aipi23 + A1P222 + P12P112 - 2pnpi22 - 2pi22A 2 + 4p223Ao + P22P111 

+ 2A 2 pll3 - 4pi33A + P112A3 

(—21) = ~pnpn3 - A1P123 - 2pi 33 A + 2P223A0 + P13P111 

( — 23) = -P11P112 - P122A1 - 4A pi23 + 2A0P222 + A1P113 + P12P111 

Proof. The simplest way to construct these relations is through cross multiplication of the 4-index relations. 
For example, substituting using equations (|4T|) and (|42|) into 



_d_ 

du 2 



(P3333) 



_d_ 

du?, 



(P2333) 







gives the first relation in the set above. 



□ 



The bilinear relations in the (3, 4)-case were considered in [13] but again, a complete set was not available 
until now. 

Theorem 4.5. Every bilinear relation associated with the cyclic (3,4)-cwi;e may be given as a linear com- 
bination of the 21 below. 



(-6 

(-7 
(-8 

(-9 

(-10 
(-11 
(-11 
(-12 
(-12 
(-13 

(-14 
(-14 

(-15 

(-15 

(-16 

(-17 



= p222 + A3P333 + 2p23p333 - 2p 33 p233 

= 2pi 33 + p22p333 + P23P233 ~ 2p 33 p 2 23 

= — 4pl23 + 4p22p233 ~ 2p23p223 ~ 2p33p222 

= 4pi22 - A 3 p222 - A 3 2 p333 + 4A 2 p333 ~ 2p 2 3p222 + 2p 2 2p223 

+ 8pi 3 p333 - 8p 33 pi33 
= -pl2p333 + P23P133 + 2p 33 pi23 ~ 2pi 3 p 2 33 
= 4pu3 - 3pi2p233 - P22P133 + 2p23pl23 + 2p33pi22 
= pll3 - 3pi 3 p223 + 2p22pl33 + 2p23pl23 ~ P33P122 
= A 3 (pi 3 p333 - P33P133) + A1P333 + pi2p223 - P23P122 

= P13P222 - Pll2 - 2A 3 (pi3p333 - P33P133) - 2A1P333 ~ 2p 2 2pl23 + P23pl22 
= 3As(pi3p233 — P23P133) — 2A2pl33 + Aip233 + 3pi2p222 

- 3p22pl22 + 4pnp333 + 4pi 3 pi 33 - 8p 33 pn3 
= -4pi2pi33 + 4pi 3 pi23 - 2p23pH3 + 2p 33 pn2 

= 5A3PH3 - 3A 2 pl23 + |Aip223 - 3pnp233 - Pl2pl33 - 2pi3pi23 

+ 4p23pll3 + 2p33pH2 
= -pllp223 - f P12P123 - §Pl3pl22 + |p22pll3 + |p23pll2 + §Aip 2 22 ~ Pill 

+ 2A3PH2 — §A2(pl3p333 — P33P133) — A2P122 — ^A3Aip333 — 3A0P333 
= — |pl2pl23 + §Pl3pl22 - 3P22PH3 + 5P23P112 

— §A2(pl3p333 — P33P133) — |Aip222 — ^A 3 Aip 33 3 — |Aop333 
= -A 3 (piip333 + pl3pl33 - 2p33pii 3 ) - 2A 2 (pl3p233 ~ P23P133) 

+ 3Aipi33 - 4A p233 - Pllp222 - P12P122 + 2p22pll2 
= f A2P113 - f A1P123 + 2A p223 - |pllPl33 + §P13P113 + §P33P111 
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( — 18) = -A 3 A p333 + A p222 + f A 2 pii2 - |pllpl23 + 3P12PII3 

+ |pl3pll2 + |p23plll - |AlPl22 ~ Ai(pi3p 33 3 - P33P133) 

(—18) = -A3A0P333 - 2A p222 - 3A2P112 - 3P12P113 + 3P11P123 

+ 3P13PU2 - |P23P111 + |Aipi22 ~ Al(pl3p333 - P33P133) 
( — 19) = §A p 13 3 - |A 2 2 pi33 - 2Ai (pi 3 p233 ~ P23P133) + Pl3pl33 ~ 2p 33 pll3) 

+ IA2A1P233 + IA3A1P133 — 2A 3 Aop233 — 3PHP122 + 3pl2pll2 + |p22plll 
( — 21) = 5Ai 2 p 333 - 2A 2 A p333 - 3A pl22 - 4A (pl3p333 ~ P33P133) 

+ 5A1P112 + Pl3plll - PllPll3 

(—22) = -5A2A1P133 - P11P112 + P12P111 - 4A (pi3p233 - P23P133) 

- 3-Al(pllp 3 33 + P13P133 - 2p 3 3pll3) + |Al 2 p233 + 3A 3 Aopi 3 3 - 2A 2 A p233 

Proof. Similarly to Theorem 14. 4[ these can be derived through cross-multiplication of the 4-index relations. 
However, the existence of Q1333 in the 4-indcx relations means that more care has to be taken in the choice 
of cross products. In higher genus trigonal cases, (or in the case where n > 3) the inclusion of further 
Q-functions in the basis makes this method increasingly tricky. An alternative method to systematically find 
bilinear relations has been developed and is discussed in [T7] . 

□ 

5 Addition Formulae 

Here we discuss the addition formulae satisfied by the Abelian functions and present some new formulae 
associated with automorphisms of the curves. We start by considering the formulae which generalise equation 
© from the elliptic case. Such a formula will exist for every (n, s)-curve as demonstrated by the following 
theorem. 

Theorem 5.1. Given an (n, s)-curve, the associated functions satisfy a two-term two-variable addition for- 
mula of the form 

a(u + v)cr(u — v) ^-^ , . . „ , . , „ 

1 hfU = E CiMv)Bi(v) 45 
a[uYa(yY 

where the functions Ai(u), Bi{v) belong to the basis for T(J, O(20[ 0_1 ])) and the ci are constants. Further, 
the polynomial on the right hand side will either be symmetric or anti- symmetric with respect to the change 
of variables (u,v) <— > (v,u) when the a-function is odd or even respectively. 

Proof. Denote the left hand side of equation (|4"5"]) by LHS(m, v). Clearly this has poles of order at most two 
and so we just need to prove that it is Abelian. We let £ be a point in the lattice and use the quasi-periodicity 
condition flU} with * = (?]'£' + rf£") T to check that 

t t / „ \ tr(u + £ + v)a(u + £ - v) x(£)e* lu+v+ ^o-(u + v)x(£)e^ u -- v+ ^cr(u-v) 

LHbm + £,v) = ; = — — — 

a(u + £) 2 (T(v) 2 X (£) 2 e lu+ * ] o-(u) 2 o-(v) 2 

= g(M I "l a{u - v) e nu + vH^-*H- 2 »-t\ = lhs(m, v). 

o-(u) z a(v)' ! 

Then using the parity property (|23|) with k = l/24(?i 2 — l)(s 2 — 1) we can check that 

a(u + v + £)cr(u - v - £) a(v + £ + u)(-l) k cr(v +£-u) 



LRS(u,v + £) = 



a(u) 2 a(v + t) 2 a{v + £) 2 a(u) 

l) k LHS(v + £,u) = (-l) k LRS(v 7 u) = (-1^ + 



a(v) 2 a(u)'^ 



= { _ lf Mu + v)o-(u-v) 
a[uYa{vy 



Hence the left hand side is Abelian with respect to both u and v and so the right hand side must be expressed 
using the basis elements for T(J, 0(20l 9-1 l)) as described. The further symmetry property of the right hand 
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side can be concluded by simply applying the symmetry property of o(u) to LHS(u,m). 

□ 

The coefficients in the right hand side of the formulae can be explicitly determined using the cr-expansion. 
(See [TH] for details of such calculations). In [T3] the authors showed that the functions associated with the 
cyclic (3,4)-curve satisfy 

a(u + v)a(u — v) , . , _ . , , , . , . . , . , . 

, ^2 I \2 = "Pi! u ) - 3^1333 (W)P33 N + Pl2{V)p23{u) + pl 3 « P22 «) 

+ pll(v) + \Q\ 333 {u)p 33 (v) - f»12(«)p23(«) ~ Pl3(")p22(u). (46) 

Note that the (3,4) a- function is odd and hence the addition formula here is anti-symmetric in (u,v). The 
corresponding formula for the functions associated to the cyclic (2, 7)-curve is 

a(u + v)a(u — v) 

= Am) - p n (v )p33 W - p22{u)p 13 {v) + pi2 f )p23 « + 2pi 3 (u)pl 3 («) 

o\uy<j(yy 

+ A(v) - pn(u)p 33 (v) - p22{v)p 13 {u) + p 12 {u)p 23 {v), (47) 

as first established in [SJ. This time the a- function is even and hence the formula is symmetric in (u, v). 

In some cases there are more addition formulae associated with the functions, resulting from automor- 
phisms of the curve equation. Such addition formulae were the topic of [15j which gave a thorough treatment 
of the genus one and two cases. We will present two new genus three addition formulae associated with the 
cyclic (3, 4)-curve. The first of these is related to the automorphism on the curve (|12|) given by the operator 

/ 27TI 

[C] : (x, y) i-> (x, (y), where ( = exp f — 

So C is a cube root of unity and [£] an operator which multiplies y by the root leaving the curve unchanged. 
We extend this notation to define the sequence of operators and automorphisms, 

[C 3 ] : (x,y) h-> (x,( j y), forjeZ. 

We can check using the basis of differentials (|15|) that these operators act on the variables u as follows. 

[C> = (CV,CV,C 2 ^3). (48) 

The action of such operators on the lattice A is stable, (it moves the points around but does not change 
the overall lattice) . This can be checked by considering the effect on the individual elements of the period 
matrices, (see [Mj for more details). We can now derive the following result for the a- function which follows 
Lemma 4.2.5 in [2^] . 

Lemma 5.2. The a -function associated to the cyclic (3,4)-curue satisfies 

o-([C j ]u) = C j a(u). (49) 
Proof. Consider the quasi-periodicity of a([( J ]u). If £ is a point on the lattice then 

a([C j ](u + £))=a([^]u+[C j }£). 
Since the lattice is stable under the action we know that [C, 3 ]l is also on the lattice. Hence by equation (|2~4"|) 

a([C j ](u + £)) = x([C^)<x([C>) exp [l([C 3 ] (« + I), [C]*)" 
In j24j the author shows that for an automorphism of a cyclic curve we have 

L([()u,v) = L{u, [C» 

and hence 

L([C>, K 3 ]v) = L(u, [C J ][C» = L(u,v). 
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Therefore, we have 

a([C](u + £)) = X ([C j ]£HK j }u) exp 
We now consider the quotient 



L[u + -,t 



<j(mu + e)) = x([C J ']*M[C J » = ± ^([C J » 

a((u + £)) cr(u) <r(u) 

since xCQ = ±1- So we see that the function 

a{u) 

is bounded and entire (since the zero sets coincide). Hence, by Liouville's theorem, the function is a constant. 
Using the Schur-Weierstrass polynomial ([25]) we see that this constant is C, 3 . 

□ 

We can now derive the addition formula associated with these automorphisms. Note that this is a more 
general version of the formula presented in Theorem 10.1 of [13) . 

Theorem 5.3. The functions associated to the cyclic (3, 4)-curve satisfy 

a(u + v + w)a{u + [Qv + {( 2 }w)q-(u + {( 2 }v + [Qw) 
a(u) 3 <j(v) 3 o-(w) 3 

= f(u, v, w) + f(u, w, v) + f(v, u, w) + f(v, w, u) + f(w, u, v) + f(w, v, u) 

where 

f(u, v, w) = P 3Q + P 27 + P 24 + P 21 + P 18 + P 15 + Pi2 +P 9 + P 6 + P 3 + P Q («, v, w) 
and the polynomials Pk(u,v,w) are as presented in Appendix[H\ 

Proof. Denote the left hand side of the formula by LHS(it, v, w). Using Lemma 15721 and the parity property 
of the (7-function we first check that LHS(it, v, w) is symmetric under all permutations of (u,v,w). Next 
consider the affect of u i— > u + I. 

a(u + t + v + w) a(u + t + [(]v + [£ 2 }w)<j{u + I + [( 2 }v + [Qw) 
LHS(m + I, v, w) = 



LHS(u, v, w 



a(u + £) 3 a(v) 3 a(w) 3 
X(^)e* [u+u+w+£/2l x(^)e* [ " +[cl " +[c2lw+i!/2l x(^)e* [ '" +[c2l ' u+[C1 ' ll ' +i ' /21 




x (£)3 e 3*[u+^/2] 

= LHS(«, v, t „) e *-[ 1 +K]+[C 2 ]] e *-[ 1 +[C]+[C 2 ]] 
However, from (|48|) we see that 

/ «i(i + C + C 2 ) 
^(i + [C] + [C 2 ]) - v 2 (i + c + ( 2 ) 

V v a (i+e + o 

and so LHS(it, v, w) is Abelian with respect to u. Further, since it is symmetric in (u, v, w) we can conclude 
that it is Abelian in v and w as well. Hence it may be expressed as 

LE.S(u,v,w) = J2aAi(u)Bi(v)Ci(w) 

i 

such that the c, arc constants and the functions Ai(u), Bi(v),Ci(w) belong to the basis for T(J, 0(39[ 2 1)), 
presented earlier in equation (|38p . To determine the constants Ci we use the tr-expansion. The computations 
involved can be heavy and so it is essential that we take into account all the available simplifications. We have 
already noted that LHS(u, v, w) is symmetric under all permutations of (it, d, w) and we reduce the number 
of independent a by applying this property to the sum. We can also check using the parity property of 
cr(u) that LHS(u, v, w) is even under (u, v, w) \-> [—!](«, v, w). Since we know the parity of the p-functions 
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matches that of the number of their indices, we can check the parity of all basis functions and hence only 
include suitable combinations. The biggest computational simplification come from the knowledge that 
LHS(it, v, w) has total weight —30 which will drastically reduces the number of possible terms. To further 
ease the time and memory constraints we implement code in Maple to efficiently expand the products of 
series so that only the relevant terms are considered. 

We find that LHS(tt, v, w) is given as stated in the theorem. For simplicity we group together the terms 
with common weight ratios into the polynomials Pk(u,v,w) which contain the terms with weight — fc in 
the Abelian functions and weight (—30 — k) in the curve parameters. These polynomials are presented in 
Appendix |E] and are made available in the supplementary material. 

□ 

If we were to try and derive the corresponding addition formula for the (2, 7)-case then we would be 
led to consider a curve automorphism [£] where the constant £ = — 1 instead. The automorphism addition 
formula would then coincide with the standard addition formula (|47j) . Hence there is no corresponding 
addition formula to Theorem 15.31 in the (2, 7)-case, or rather it is the same as the tradition addition formula 
in equation (j4"T)) . 

The final addition formula presented in this paper is satisfied by the functions associated to the reduction 
of the (3, 4)-curve given by y 3 = x A + Ao- This has a family of automorphisms, 

[i J ] : (x,y) i y ((-i) J x,y) 

where i is the complex variable and j S Z. The functions then satisfy the following formula. 

Theorem 5.4. The functions associated to the restricted (2>,4)-curve, y 3 = x A + Ao satisfy 

a(u + v)cr(u + [i]v)a(u + [i 2 ]v)a(u + [i 3 ]v) 



= f(u,v) - f(v,u) 



where 



f(u,v) = -|pll 2 (v)pl23N - gPllllH + J E pll2(u)d 2 p lll] (v) - \F{u)p 2 2{v) 

+ 57P1113(M)P3333('W) ~ jgd 2 p lll] (u)d 2 p ll3] (v) - ^d 2 d 2 Ql333(u)p2222{v) 

+ \d2p [33] {u)p2Z3{v) - §PU13(m)P22(u) + fp^ 1 («)Ql333(«) ~ § £>X133(«)Ql333(«) 

+ ^9 2 a 2 Ql333(u)pl333(«) - JsP3333{v)F(u) + | p [22] («)pl333 («) 

- f<9 2 p[13](»Pl23(» - ^Pi133(m)P2222(«) + [\p2222{v) + jg p 3 333 («) P22 («)] A . 

Proof. Using a similar approach to Lemma [5.21 we see that cr([i]it) = —ia(u). We can then check that the 
left hand side is anti-symmetric under (u, v) i— > (v, u) and is Abelian. We can hence express it as a quadratic 
polynomial in the elements of the basis for r(j, 0(4Ol 2 !)), derived earlier in Theorem (|3.7[) . 

We use the cr-function to determine the coefficients of each term. Due to the restriction on the curve 
parameters, determining the coefficients is computationally easy in comparison to Theorem 15.31 (However 
the construction of the basis in Theorem p.7[) required much effort.) 

□ 

There is a similar formula associated to the restricted (2, 7)-curve y 2 = x 7 + Xq which has automorphisms 
[P] : (x, y) \-¥ (P x, y). Here l is a seventh root of unity and we consider, 

nL, 

a(u) 7 a(v) 7 

To evaluate this we will require a basis for the 7-pole vector space which has dimension 7 3 = 343. 
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A Bases for the functions associated with the hyperelliptic curve 
of genus four 

In this Appendix we consider the general (2,9)-curve, 



fl 3 X 



/j, 5 x 
„9 



H 7 x 

„8 



w)y 

„7 



X + jJL 2 X + fHX + [1§X + fi S X + filQX + [l\ 2 X +/J14I + fliQX + /lis, 



which has genus four. We construct the bases for standard Abelian functions associated with this curve, 
following the approach in Section [5J 



Theorem A.l. The basis for C(29[ 3 1)) is given by 



© CI 

© Cp 23 

© CA 2 



© Cp 2 4 

© CA 3 



Cpi2 © Cpi 3 
Cp 33 © Cp 34 

CA 4 © CA 5 , 



H4 



^P22 

CAi 



(50) 



where 



Ai 
A 2 
A 3 
A 4 
A 5 



P34pl4 



P24 



P33p24 + P44P13 ~ P22p44, 



P34P23 

P34P13 + P24P14 - P33P14 ~ Pl2p44, 
-p44pll + P14 - P23P14 + P13P24, 
-2p 3 4p n + 2p 13 p 14 - 2p 22 pl4 + 2pi 2 p24, 
-P12P23 + P22P13 - P13 + Pl2pl4 ~ Pllp24 



Pllp33- 



Proof. We follow the proof of Theorem 13.31 This time the dimension is 2 9 — 2 4 — 16 and we have ten 
2-index p-functions. By testing arbitrary sums of quadratic terms in the 2-index p-functions we find that we 
can identify combinations that have poles of order only two. We find five linearly independent combinations 
which can fill the missing basis entries. 

□ 

Note that we could have alternatively used Q-functions as discussed in Section 13.11 However, The A- 
functions are advantageous since it allows the theory to be completely realised in terms of 2 and 3-indcx 
p-functions. As discussed in Section 21 this appears to be a feature unique to the hyperelliptic cases. Note 
that while the T, J- and ^/-functions introduced in Section [3] had reduced poles structures in general, these 
A-function have poles of order two only in the (2,9)-case. 

Theorem A.2. The basis for T(j, C(39[ 3 ])) is given by 



where 



(37J © Cpin © 


Cpn2 9 


) Cpii 3 


© Cpn 4 


© Cpi22 


© Cpi23 © 


Cpi24 9 


) Cpl33 


© Cp 434 


© Cp i44 


© Cp 222 © 


Cp223 S 


) Cp224 


© Cp233 


© Cp 23 4 


© Cp 2 44 © 


Cp333 9 


3 Cp334 


© Cp 34 4 


© Cp444 


© CdiAi © 


Cd 2 Ai 9 


) O3A1 


© O4A1 


© C%A 2 


© C<9 4 A 2 © 


OiAa 9 


5 C5 2 A 3 


© Cd 3 A 3 


© G9 4 A 3 


© C<9iA 4 © 


C9 2 A 4 9 


3 Cd 3 A 4 


© C9 4 A 4 


© C9iA 5 


© C<9 2 A 5 © 


C9 3 A 5 9 


3 C5 4 A 5 


© C7lll333 


© C7Ill334 


© CT111344 © 


C7ii2334 9 


3 C7l22244 


© C7i 2 2334 


© C7i22344 


© CT122444 © 


C7i23333 9 


3 C7i23334 


© C7i23444 


© C7i 24444 


© CT133333 © 


C7i44444 9 


3 C7222344 


© C7222444 


© C7223444 


© C7^24444 © 


C7233333 9 


3 C7233334 


© C7233444 


© C7333334 


© C7333444 © 


C7334444 9 


3 C 7344444 


© CUi 


© cu 2 . 


2p233p2223 — p222p2333 


3 

— 2 p223p2233 


— P33P22223 


+ p22p22333 + 


5P333P2222 


+ 6p 22 p33p223 + 6p 23 p33p222 + 3p| 3 


p223 - 12p22p23p233 ~ 3pf; 


>P333i 


2p234p2223 — p222p2334 


— |p224p2233 


— P223P2234 


— P34P22223 + 


P22p22334 



(51) 



;p334p2222 + 6p 2 2p34p223 + 6p 2 3p34p222 + 3p2 3 p224 " 12p22p23p234 ~ 3p| 2 p 334 . 
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Proof. This time the dimension is 3 9 = 3 4 = 81 and we find that we can identify 54 functions using the basis 
(|3"2"j) and its derivatives. We then proceed to add 36 of the T- functions defined in equation (|36|) . To find 
the final two functions we used a new class formed from a sum of p-functions in which each term has seven 
indices. 

□ 

In the (2,9)-case tr(it) has weight 10 and hence the minimal weight for functions in this basis is 30, in 
accordance with Lemma 13.61 This is achieved by the function 7iii333, although this is not a unique choice. 
For example, it could be replaced by 7222222 ■ Both functions will reduce to a constant over a(u) 3 . 

B The 4-index relations associated to genus three curves 

The complete set of 4-indcx relations associated to the (2, 7)-curve is given below. 

(-4) p3333 = 4p 23 + 4p 33 A 6 + 2A 5 + 6P33 

( — 6) P2333 = 6pi3 - 2p 22 + 4p 23 A 6 + 6p23p33 

(—8) P2233 = ~2pi2 + 4pi 3 A 6 + 2p 23 A 5 + 2p 22 p33 + 4p 23 
(—8) P1333 = -2pi2 + 4pi 3 A 6 + 6pi3p 33 
( — 10) P2223 = -6pn + 4pi 3 A 5 + 4p 23 A 4 - 2p 33 A 3 - 4A 2 + 6p 2 2p23 

( — 10) pi 233 = 2pi 3 A 5 + 2pi2p 3 3 + 4pi 3 p 23 

( — 12) pi 223 = 2pi2p23 + 4pi 3 p22 + 4p i3 A 4 - 2Ai - 2 pl 3 + 2pnp 33 
( — 12) P2222 = 6p22 - 6A1 + 12pi2p23 + 12p? 3 - 12pi 3 p22 - 12pnp 33 + 2A 5 A 3 - 8A 6 A 2 
+ 4P22A4 - 12p 33 A2 + 4p i2 A5 + 4p 23 A 3 - 12piiAg 

( — 12) pll33 = 2pl2p23 + 6pi 3 - 2pi 3 p 22 

(—14) pn2 3 = 2piip 23 + 4pi2pi 3 + 2pi 3 A 3 - 4A 

(—14) P1222 = 6P12P22 - 2pnA 5 + 4pi 3 A 3 - 6p 33 Ai + 4pi 2 A 4 - 8A - 4A 6 Ai 

(—16) pm 3 = -2p 23 Ai + 4p 33 A + 4p 13 A 2 + 6p u pi 3 

(—16) P1122 = -2P23A1 + 2pi 2 A 3 - 8p 33 A + 4p 43 A 2 + 2pnp 22 + 4pJ 2 - 8A 6 A 

(—18) pni2 = 6pnpi2 - 4A 5 A - 2p 22 Ai + 6pi 3 Ai - 8p 23 A + 4pi 2 A 2 

(-20) p im = 6p n + 2A 3 Ai - 8A 4 A + 4p u A 2 + 4p 12 Ai - 12p 22 A + 16pi 3 A 

The complete set of 4-index relations associated to the (3, 4)-curve is given below. 

(-4) p3333 = -3 p22 + 6p 33 

( — 5) p2333 = 3p33A 3 + 6p 23 p 33 

(—6) p 22 33 = 4pi 3 + 3p 23 A 3 + 2A 2 + 2p 22 p 33 + 4p 23 

( — 7) p 2 223 = 3p 22 A 3 + 6p 22 p23 

(—8) P2222 = -3P33A3 2 + 12p 33 A 2 - 4Qi 333 + 6p 22 



(-8) 


P1333 


= Ql333 + 6pi3p 33 


(-9) 


P1233 


= 3pi 3 A 3 + A 4 + 2pi 2 p 33 + 4pi 3 p 23 


-10) 


P1223 


= -2pn + 3p 42 A 3 + 4p 42 p 23 + 2pi 3 p 22 


-11) 


P1222 


= -Ql333 A3 + 6p 33 Ai + 6pi 2 p22 


-12) 


P1133 


= 2pi 3 A 2 - P23A1 + 2pnp 33 + 4p^ 3 


-13) 


P1123 


= 2pi 2 A 2 - P22A1 + 2piip 23 + 4pi2pi3 


-14) 


P1122 


= -IQ1333A2 + P33A3A1 + 8p 33 A + 2pnp22 + 4p^ 2 


-16) 


P1113 


= 3p i2 Ai - 6p 22 A + 6pnpi3 


-17) 


P1112 


= -Q1333A1 + 6p3 3 A 3 Ao + 6pnpi 2 


-20) 


P1111 


= -4Q1333A0 - 3p 33 Ai 2 + 12p 33 A 2 A + 6p n 
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C Quadratic 3-index relations associated with the cyclic hyperel- 
liptic curve of genus three 

This appendix contains the complete set of quadratic 3-index relations associated to the cyclic (2, 7)-curve. 
Note that A is a function quadratic in the 2-index p- functions, as defined in equation (|33[) . It was used 
in the basis of fundamental Abelian functions (|32j) . and occurs here only linearly, or multiplied by a single 
2-index p-function. Hence each relation can be rewritten as a polynomial in 2-index p-functions of degree 
three. The relations are presented in weight order as indicated by the number in brackets. 



(-6) 
(-8) 

(-10) 
(-10) 

(-10) 
(-12) 

(-12) 



(-12) 
(-12) 

(-14) 
(-14) 



4 P 33 + 4A 4 + 4A 5 p 33 + 4A 6 p§ 3 - 4pi 3 + 4p 2 2 + 4p 33 p 2 3 
P233P333 = 4p23p33 + 2 A 3 + 2A 5 p 23 + 4A 6 p 33 p23 + 2pi2 + 4p 33 pi 3 



2 

P333 



2 

P233 
P223P333 



P133P333 
P223P233 

P222P333 



P133P233 
P123P333 



P223 



P133P223 



= 2A 



■2A 5 p|3-4A 2 p33 + 2Af 



-16p? 



:2A 

: 2pnp33 



" 2p 3 3p22 + 2p| 3 

P23P13 - 4p23p22 + 4pn + 4A 2 + 4p| 3 p 33 + 4A 6 p§ 3 

"4p 3 3pl2 - 2p23pl3 + 4p23p22 + 4A 6 p33pl3 + 2A 5 p 33 p2 3 + 4A4p 23 

- 2pu + 2A 5 pi 3 - 2A 3 p 33 + 2p22p| 3 + 2p^ 3 p 33 

-2p33pi2 + 2p2 3 pi3 + 4A 6 p 33 pi 3 - 2pn + 2A 5 pi 3 + 4pi 3 p| 3 

1 + 4A 6 p 23 pl 3 + 2p22p23p33 + 2p^ 3 - 2pi 3 p 22 + 4p^ 3 - 4pnp 33 
f 2A 3 p23 

13 _ *¥<2a. T 20pi 3 p22 - 2pi2p 23 + 6p 22 p23p33 + 4piip 33 - 2p| 3 
2A5A3 - 8A 6 p23pl3 - 2A5P23 - 2A 5 pi2 - 4A 4 p22 + I6A4P13 - 4A 3 p 3 
8A4p 33 p 23 - 2A 5 p 33 p22 - 4A 6 p 33 pl2 + 8A 5 p 33 pi 3 + 8A 6 p23p22 

8A 6 A 4 p23 - 4A 6 A 3 p 33 - 4A 5 p| 3 - 2A 3 p 23 

4A 6 p 23 pl 3 + 4pi3p23p33 - 2pi 3 p22 ~ 2pi 2 p23 + 4p 
- 4 Pl3 + 4pl3p22 + 4A 4 pi 3 + 2pi2pi 3 
+ 2pi 3 p 23 p 33 + 2A 5 p33pi 3 

4p 33 A - 4p 23 pn - 4pi 3 pi2 + 4A 
+ 4A 5 p 23 pi3 + 4A 2 p 33 + 4A '- 



>? 3 



-4A 



>6Pi 3 



4A 3 p33p23 
<4p| 3 - 4Aip 33 + 4p22pl 3 
= 4A - 2p 23 pii - 4pi 3 pi2 + 4A 6 p^ 3 + 2A 5 p 2 3pl3 - 2Aip 33 

+ 2pi 3 pi 3 + 2pi 3 p 33 p 22 
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14) 
14) 



14) 
14) 



14) 
16) 



16) 



16) 

16) 
16) 



16) 
16) 



18) 



18) 



Pl23p233 
Pl22p333 



Pl33 = 4 P13P33 - 4pi 3 pl2 + 4A + 4A 6 p 2 3 
P222P233 = -8p 33 A - 4p23pn + 4A 6 A 3 p 2 3 - 4A 6 A - 4A 6 p 2 3 - 8A 6 A 2 p 3 3 

+ 2A 5 p23p22 + 8pi 3 pi2 - 2p 22 pl2 ~ 4A 2 A 5 + 4A 3 p 33 p23 + 2p| 2 P33 

+ 2p22pi 3 - 8A 2 p 2 3 + 4A 6 pi 3 p 22 - 4A 6 pnp 33 - 4Aip 33 - 4A 5 pn 
+ 8A 3 pi 3 - 2A 3 p22 - 4A 2 p2 3 

2pl3pl2 + 2A 5 p 23 pi 3 - 2Aip 33 + 2pi 3 p2 3 + 2pi2p 33 p23 + 2A 3 pi 3 

4p 33 A - 2A 5 p 33 pi2 - 8A 6 p 2 3 + 8A 6 A 4 pi 3 + 8A 4 p 33 pi 3 + 6pi 3 pi 2 

- 4A 5 p 23 pi 3 - 2A 5 2 pi 3 - 4A 4 pi2 + 2A 5 pn + 2A 3 p i3 + 8A 6 pi 3 p 2 2 

+ 4A 6 pnp 33 - 2pi 3 p| 3 + 2pi 3 p 33 p 22 + 4pi2p 33 p 23 - 4p 22 pl2 + 4p 23 pn 
P113P333 = -2p33A + 2p 23 pll + 2pi 3 pi 2 + 2p 2 3 p 33 + 2pup| 3 + 2A 3 pi 3 
P222P223 = -4A 6 p23pll — 4p 2 3A + 4A 5 pi 2 p23 ~ 4A 5 pnp 33 + 4A 4 p 23 p22 

- 2A 3 p 33 p22 - 4A 6 Aip 33 + 2A 5 A 3 p 23 - 4A 2 p 33 p2 3 - 4A 2 A 5 p 33 - 8A 4 pn 
+ 2A 3 p2 3 + 4A 3 pi2 + 8A 2 pi 3 - 8A2P22 - 2Aip 23 - 8A p 33 

+ 4pi 2 p23 + 4A 5 pf 3 - 6P22P11 + 2p 2 2 + 2A 3 2 - 8A 2 A 4 - 2A 5 Ai 
P133P222 = -4A 6 p23pn - 4p23A + 2A 5 pi2p23 - 4A 6 Aip 33 - 2pi 2 p| 3 

- 2Aip 23 - 8A0P33 + 4A 5 p 2 3 - 8pi 3 pn + 2p22pn 

- 2A 5 Ai + 4pi 3 p 23 p22 + 2pi2p 33 p22 - 4Aip| 3 

2p 23 A + 2pi2p2 3 - 4A p 33 - 2pi 3 pn + 2A 5 p 2 3 - 2A 3 p 33 p 43 
+ 4A 4 p 23 pi 3 + 2pi 3 p 23 p22 + 2Aip^ 3 
-4A p 33 - 2pi 3 pn + 2p 2 3 p 23 + 2pi2pi 3 p 33 + 2A 5 p 2 3 
2A 5 pi2p2 3 - 2A 5 pnp 33 - 4A 6 Aip 33 + 2p 42 p2 3 - 2A 3 p 42 + 4A 2 pi 3 

- 4Aip 23 - 2A 5 p 2 3 + 4pi 3 pn - 2p 2 2 - 2A 5 A 4 + 4A 6 pi 3 pi 2 
+ 4A 6 A 3 pi 3 + 4A 3 p 33 pi 3 + 2pi2p 33 p22 + 2A 5 A - 4Aip| 3 



P123P223 

P123P133 
Pl22p233 



P113P233 
Pll2p333 



2 

P222 



-2p 23 A + 4A 2 pi3 - 2A1P23 + 4pi 3 pii + 2p 13 p 23 + 2piip 23 p33 
4A 6 p2 3 pii + 2p 23 A - 2A 5 pnp 33 - 2A 3 pi2 - 2p 2 3 p 23 - 2A 5 p 2 3 
+ 4pi 3 pn - 2p22pn - 2p 2 2 + 4A 6 pi 3 pi2 + 4A 6 A 3 pi 3 
+ 4A 3 p 33 pi 3 + 4pi2pi 3 p 33 + 2pnp 23 p 33 + 2A 5 A 
-8A 3 p 23 pi 3 + 8p 23 A 5 A 2 - 8A 5 A 3 pi 3 + 4p 23 p22A 3 - 16A 6 A 4 pn 
- 16A 6 A 2 p22 + 4A 5 A 3 p 2 2 - 16A 2 A 4 p 33 - 8p 33 AiA 5 + 8A 3 pi 2 A 6 



16A p 33 - 4p 2 2 p 33 - 4pnp| 3 



JP13P23 - 8Aip 33 p 23 



4pl 2 



18) 



+ 16A 6 pi 3 pn + 16A 6 A 2 pi 3 + 4piip 22 p 33 - 4p 2 3 p 22 + 16Aip 43 

- 16AiA 4 - 16A 6 A 4 A 2 - 16p 2 2pnA 6 - 16p 2 2A 2 p 33 + 4A 5 p 22 pi2 
+ 4A 6 A 3 2 + 4A 5 2 A 2 + 4p| 3 A 2 + 4p 2 2 A 4 + 4p 2 2 A 6 + 4A 5 2 pn 
+ 4A 3 2 p 33 - 16A 4 A - 8A 3 pn + 32p 43 A - 20p 22 A - 8A 5 pi 3 pi 2 

- I6A1P22 + 16A 2 p 33 pi 3 - 8pi2pn + 16A 4 p 2 3 

P123P222 = -2P33A1A5 + 4p 43 p22A 4 + 4A 4 pnp 33 - 2p 23 A 5 pii - 4Aip 2 2 

+ 8p 43 A - 2p 22 A - 2pi2pn + 4pi2pi 3 p 23 - 2Aip 33 p23 + 2pnp22p33 

+ 4A 4 p 2 3 + 8A p| 3 - 2p 2 2 p 33 - 2pnp^ 3 - 2A 3 p 33 pi2 + 2pi 3 p| 2 

- 2A 3 pn - 4A 4 A - 2p 2 3 p22 + 4Aip 43 - 4AiA 4 + 2pi 2 p2 3 p22 
P123 = 4 Pi 3 A + 4pi2pi 3 p23 + 4A 4 p 2 3 + 4A pl 3 
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P122P223 = ~4AiA 4 - 2p 33 AiA 5 + 4A 5 pi 3 pi2 + 4A 3 p23pl3 + 4pi3p22A4 - 4Aip22 

- 8A 6 A p33 - 4A 6 pi3pn - 4A 5 A - 12pi 3 A + 4p 22 A - 4pi 2 pi 3 p 2 3 

+ 2p 2 3 p 22 + 6A1P13 - 8A0P23 - 8A 4 p 2 3 - 8A p 2 3 + 2p 2 2 p 33 + 2p u p 2 . 3 

- 4A 2 p33pi3 - 2pnp 22 p33 - 2A 3 pn + 2A 5 A 3 pi3 + 4pi2p 23 p22 

P122P133 = -4pi 3 A + 2pi2pn - 4A 5 A + 2p 2 2 p 33 + 2p 2 3 p 22 

- 8A 6 A p33 + 2A 5 pi 3 pi2 - 4A 6 pi3pn + 2Aipi 3 - 8A0P23 - 8A p| 3 

Pll3p223 = -4pi 3 A + 2pnp| 3 + 2A 3 p 23 pl3 + 2Aip 33 p23 ~ 4A 2 p 3 3pl3 + 2p 2 gp22 

+ 2Aipi 3 - 4A p23 
Pll3pl33 = -2pl3A + 2piipi3p 33 + 2pf 3 + 2Aipi 3 - 4A p23 

P112P233 = 8pi 3 A - 2p 22 A - 4pi2pn + 4pi2pl3p23 - 4Aip 33 p23 + 8A 2 p33pl3 

+ 8A 6 pi3pn + 8A 6 A 2 pi3 + 2pnp22p33 - 4A 6 Aip23 - 4A 2 pi2 

- 2 p? 3 p22 - 4Aipi3 + 2Aip22 

Plllp333 = 2A 5 pi 3 pi2 + 2A 3 p23pl3 + 6pnpi3p33 - 4A 4 piip3 3 + 2p2 3 A5pii 

- 4pi2pl 3 p23 - 2pnp22p33 + 2p 2 3 p 2 2 - 4A 4 p 2 3 + 2p 2 2 p 33 + 2pnp| 3 
+ 2A 5 A 3 pi3 + 4A 4 A - 10pi 3 A + 4p 22 A - 2p? 3 + 2A3P33P12 

P122P222 = -8A 6 A 4 Ai - 4A2P33P12 - 6A1P33P22 + 4A3P13P22 - I6A0P22 + 2Aip 23 
+ 4A 5 Aip23 + 8A 6 Aipi 3 - 8A5A0P33 + 4A 4 p 22 pi2 - 8AiA 4 p 33 

- 4A 6 pi2pn - 8A 6 Aip 22 - 4A 6 A 3 pn + 4pi 2 p 22 + 2A 5 p 2 2 - 4p 42 A 

- 16A 4 A - 2Aip 42 - 4A 3 A - 4A 2 pn + 2A 5 A 3 pi 2 + 16A pi 3 + 2A 5 2 Ai 

- 2A 5 p22pn - 2A3A1 

P122P123 = 2pi 2 p23 - 8A 4 A + 2p i2 A + 2pfj + 2pi 2 p 2 2pi3 + 4A 4 pi 3 pi 2 

- 2A 5 pi 3 pii - 2Aip 33 pi3 - 4A 5 A p33 - 8A0P22 + 2A 3 p 2 3 + 8A pi3 
P113P222 = 2A 3 pnp33 + 8A0P33P23 - 2pf 2 p23 - 2A 3 Ai - 4p i2 A - 4pf x 

- 8Aip 33 pi3 - 4A 2 p33pl2 + 2Aip 33 p22 + 2A 3 pi 3 p 2 2 + 2A 3 p 2 3 - 2A 4 pi2 
+ 4pi2p22pl3 + 2pnp22p23 - 2A 3 A - 4A 2 pll 

P113P123 = 4A P33P23 + 2p 2 3 pi 2 + 2piipi 3 p 23 - 2Aip 33 pi 3 + 2A 3 p 2 3 
P112P223 = -8A p33p23 + 2pf 2 p23 ~ 2A 3 A 4 - 4pfj + 2pnp22p23 + 4A 5 pi3pn 

+ 4Aip 33 pi3 + 4A 2 A 5 pi3 - 8A 6 A p23 + 4A 2 p23pl3 + 2A 3 pi 3 p22 - 4A 2 pll 

+ 4A 6 Aipi3 + 4A p22 - 2Aip 2 ! 3 ~ 4A 3 p 2 3 - 4Aipi 2 - 8A pi 3 - 2A 5 Aip 2 3 

P112P133 = -8A p33p23 + 2pi 3 pi2 - 2p 42 A + 2pnp33pi2 + 4Aip 33 pi3 

- 8A 6 A p23 + 4A 6 Aipi3 + 4A p22 - 2Aip i2 - 8A pi3 

Plllp233 = -2A 3 piip33 - 2pi 3 pi2 + 2p 42 A + 4pnpi3p23 + 2pnp33pi2 

+ 4A 2 A 5 pi 3 + 4A 2 p23pl3 + 4A 2 p33pi2 - 2\ip 33 p22 - 2A 5 Aip23 

+ 4A 5 pi 3 pii + 4Aip 33 pi3 - 2Aip| 3 - 2A 3 p 2 3 + 2A 3 A 
P122 = -16A 6 A 4 A + 4p 2 2 p 22 + 4pnA + 4A 5 2 A - 8A 3 A - 4A 5 pi 2 pn 

- 4Aipn + 4A 2 p 2 3 - 16A 4 A p 33 - 4Aip 2 3pi3 + 8A5A0P23 

- 8A p33p22 - 8A0P12 + 4A 3 pi 3 pi2 - 4Aip 33 pi2 + 4A6P 2 ! 

+ 8A p| 3 + 16A 6 A pi3 + 4A 4 p 2 2 - 16A 6 A p 2 2 

P113P122 = 2pi3p 2 2 - 4A 3 A + 2pnpi3p 22 - 8A p33pl3 + 4A p| 3 - 2Aipn 

+ 4A 2 p 2 3 - 4Aip23pl3 + 4A p33p22 ~ 4A pl2 + 2A 3 pi 3 pi2 - 2Aip 33 pi2 
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Pll2pl23 



Plllp223 



PlllPl33 
Pll2pl22 

Pll2pll3 
Plllp222 



Pll3 = 4 PllPl3 + 4A 0p| 3 + 4A 2P? 3 - 4Aip 23 pl3 

P112P222 = 2pup22 + 2pf 2 p22 - 4A 6 A 3 Ai - 8puA - 8A 3 A - 8A 2 pnp3 3 

- 4A 3 Aip 33 - 8A 6 A 2 pn + 2A3P22P12 - 2A 5 Aip 2 2 + 8A 5 Aipi 3 

- 8A0P23 _ 4Aipn + 8Aip 23 pi 3 - 8A5A0P23 - 8A pi2 - 4Aip 33 pi2 

+ 4A 2 pi2p23 + 4A 5 pi2pn - 2Aip 23 p22 - 4A 6 Aipi 2 - &\ep\i + 4A 2 A 5 pi 2 
2pi3p? 2 - 4A 3 A + 2pnp 23 pi2 + 2A 5 Aipi 3 - 4A p| 3 

- 2Aipn + 2Aip 23 pi 3 - 4A 5 A p2 3 - 4A pi2 + 2A 3 pi 3 pi 2 
-2pi3pi2 + 4pnA + 2pnpi 3 p22 + 4pnp 23 pi2 - 4A 2 pnp 33 - 2A 3 2 pi 3 
+ 2A 3 Aip 33 + 8A 4 pi 3 pn - 2A 3 p 23 pn - 4AiA 4 p 23 + 2A 5 Aipi 3 
+ 6Aip 23 pi 3 - 4A 5 A p2 3 - 4A p33p22 - 4A 3 pi 3 pi2 + 4Aip 33 pi2 
+ 8A 2 A 4 pi 3 + 8A 2 A - 4Aip 23 p22 + 8A p 33 pi 3 + 8A2P12P23 - 4A p2 3 
: 2 PiiP33 - 2pnA + 8A0P33P13 + 2piipi 3 + 2A 5 Aipi 3 

_ 4A p2 3 + 2Aip 23 pi 3 - 4A 5 A p2 3 - 4A p 3 3p22 + 2Aip 33 pi 2 
2pnpi2p22 + 2Ai 2 - 8A 2 A - 8A 6 A 3 A + 2pf 2 - 8A 3 A p 3 3 
+ 2A5A1P12 + 4Aip 2 3 - 8A 6 A pi2 + 4A 2 pi 3 pi2 + 8A 5 A pi 3 - 4Aipnp 33 

- 2Aipi 3 p22 - 8A p 33 pi2 - 4A 6 Aipu - 8A pn - 4A 5 A p22 + 2A 3 p 2 2 
4pnpi2pi 3 + 2Ai 2 - 8A 2 A + 4A p2 3 p22 + 2AiA 

- 8A0P23P13 + 6Aip 2 3 + 4A 2 pi 3 pi2 - 2Aipnp 33 - 4Aipi 3 p 2 2 - 8A pn 
: 18Aipi 3 p 22 + 8A 3 pi 3 pn + 2Aipnp 33 + 6pnpi2p22 - I6A4A0P23 

+ 8A 4 pi2pii + 8A 3 A p 33 - 14Aip 2 3 - 4A 3 p 2 2 + 2AiA - 4A 5 p 2 ! 

- 16A p23p22 + 16A p23pl3 + 16AiA 4 pi 3 + 8A p 33 pl2 - 2A5A1P12 

- 2A3A1P23 - 4A 2 A 5 pii - 4A 2 p2 3 pn + 8A 2 A 4 pi2 - 2A 3 p22pn - 2A 3 2 p i2 

- 4AiA 4 p 2 2 + 8A2P22P12 - 4Aip 2 2 ~ 8A 2 pi 3 pi2 - 2A 5 A 3 Ai - 2pf 2 

PlllPl23 = 2p 2 !p23 + 2pnpi2pl3 - 8A p23p22 + 8A p23pl3 + 4AiA 4 pi 3 

+ 4A 3 A p 33 - 4Aip 2 3 + 2Aipnp 33 + 4Aipi 3 p 22 - 8A4A0P23 
+ 2A3pi3pn + 4A p 33 pi2 
P112 = -16A A + 4A 6 Ai 2 - 16A 6 A 2 A + 4pnp 2 2 + 8Aipi 3 pi 2 - 16A 6 A pn 

- I6A0P12P23 - 16A 2 A p33 + 4A p22 - 4Aip 22 pi2 + 4A 2 p 2 2 + 4Ai 2 p 33 
P111P122 = 2p 2 4 p22 - 4A 5 A 3 A + 2piip 2 2 - 8A 3 A p2 3 - 4Aip 23 pn + 4A 2 pi 3 pn 

+ 8A pnp 33 - 2Aipi 3 p i2 + 4AiA 4 pi2 - 2A 5 Aipn - 4A 5 A pi2 

- 8A4A0P22 - 8A pi 2 p2 3 + 8A 2 A p 33 - 8A0P22 + 4Aip 22 pi2 - 2Ai 2 p 33 
+ 24A pi 3 p22 + 16A 4 A pi 3 + 2A 3 Aipi 3 + 2A 3 pi 2 pn - 16A p 2 3 
4pnpi 3 - 4A 3 A p2 3 - 2Aip 23 pn + 4A 2 pi 3 pn + 8A pnp 33 
+ 2Aipi 3 pi2 + 2A 3 A 4 pi 3 - 4A pi2p23 + 8A 2 A p 33 - 2Ai 2 p 33 

: 2A 5 Ai 2 - 8A2A5A0 + 4p 2 x pi2 - 8A p2 3 pll + 2A 3 Aipi2 - 4A p22pl2 

- 8A 2 A p2 3 + 8A pi 3 pi2 + 4Aipi 3 pn - 4A 3 A p22 + 4A 2 pi2pn 

- 2A1P22P11 + 8A 3 A pi 3 - 8A5A0P11 + 2Aip 2 2 + 2Ai 2 p 23 
4pf x + 4Ai 2 A 4 + 4A 3 2 A - 16A 2 A 4 A + 8A 3 A pi 2 + 4A p 2 2 + 4A 2 p 2 4 

- 16A 2 A p22 + 16A pi 3 pn - I6A0P22P11 + 4Aipi2pn + 16A 2 A pi 3 
+ 4A 3 Aipn - 16A 4 A pn - 4Ai 2 p 43 + 4Ai 2 p 2 2 



P111P113 



P111P112 



Pm 
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D Quadratic 3-index relations associated with the cyclic trigonal 
curve of genus three 

This appendix contains the complete set of quadratic 3-index relations associated to the cyclic (3,4)-curve. 
Note that Q1333 was the sole Q-nmctk>n used in the basis of fundamental Abelian functions, (|37j) . It appears 
here as a linear term, or multiplied by a single 2-index p-function. The relations are presented in decreasing 
weight order as indicated by the number in brackets. 



-6) 
-7) 
-8) 
-8) 
-9) 

-9) 
-10) 

;-io) 

-10) 
-11) 

-11) 
-11) 

-12) 



-12) 
-12) 
-12) 
-13) 
-13) 
-13) 
-14) 



P333 

P233P333 
P233 
P223P333 
P223p233 

P222P333 
2 

P223 
P222P233 

P133P333 
P222P223 

P133P233 



2 

P222 



P133P223 



P123P233 



Pl22p333 



P133P222 



P123P223 



Pl22p233 



Pl33 



: 4 P! 3 + 4 Pl3 + Pis - 4 P22P33 

4p23pl 3 ~ 2 Pl2 - P22P23 + 2A 3 p§ 3 
: 4 p| 3 p33 + P22 + 4A 3 p23p33 + 4p 33 A 2 ~ |Ql333 

2 P23P33 + 2p22pl 3 - 2p^ 2 + A 3 p 2 3p33 + 4pi 3 p33 + §Ql333 

2pi 3 A 3 + 2Ai + 4p 13 p23 + 2p23p22p33 + 2p| 3 

+ P22P33A3 + 2p2 3 A 3 + 2p 23 A 2 

6p23p22p33 - 4pi2p33 ~ 2pf 3 + 4p 22 p33A 3 ~ 8pi 3 p 2 3 - p| 3 A 3 
4p2 3 p22 + 4pu - 4pi2p23 + 4pi 3 p22 ~ 4pi 2 A 3 
+ 4A 3 p22p23 + Xlpl 3 + 4A 2 p22 ~ 4A 2 p| 3 + |p33<9l333 
2p| 3 p22 + 2p| 2 p33 + 4pi2p23 + 2pi 2 A 3 + A 3 p 22 p23 

- 2A^p|j + 8A 2 p^ 3 - f P33Q1333 

P12P23 - 2pi 3 p22 + 4pi 3 p| 3 + §p33<2l333 
2p| 2 A3 - A^p23p33 + 4p2 2 p23 + 4A 2 p23p33 + 4p 3 3Ai 
|p23<3l333 

2A 3 pi 3 p33 + 2p33Ai + pi2p22 + 4p 33 p23pl3 + §p23<3l333 



4pl3A 3 



IA3Q1333 



P123P333 = 2P12P33 - 2pi 2 p 2 2 + 2p 33 p 2 3pl3 + 2A 3 pi3p 33 ~ 5p23<2l333 
4p22 + 8pnp33 - 8p? 3 - 4p 33 pi 2 A3 + 4p23pl 3 A 3 + 4pi 3 A| 

- 4p 22 p33A| + p| 3 A§ - 8pi 3 A 2 + 16p 33 A 2 p22 - 4p| 3 A 2 



- 8A - 4p22<3l333 - 4p23Ai 

2pl3p33p22 + 2pi 3 p| 3 + 2piip 33 + 2pf 3 - P33P12A3 

+ 2pi 3 A 2 + 2p2 3 pi 3 A 3 + 2A + |p22<3l333 

2pi 3 A 2 + 2pf 3 + 2p 3 3pi2A 3 + 2p 23 p33pl2 + 2pi 3 p| 3 

- 2pnp33 + 2p23pi3A 3 + 2A - |p22Ql333 
2P13P33P22 - 2pi 3 p| 3 + 4p 23 p33pl2 ~ 2pnp 33 - 6p? 3 
+ 2p 33 pi2A 3 - P23P13A3 - 2pi 3 A 2 + P23A1 + |p22<3l333 
-2p 12 p2 3 + 4p23pi 3 p22 ~ P23P12A3 + 2p 2 2pl 3 A 3 

+ 2p 33 pl2p22 + 4p 33 Ai - |p33A3<5l333 
: 2p 22 Ai + 2p22pl3A 3 - 2p§ 3 Al + 2pi2p2 3 + 2p23pl3p22 

- 2pnp23 + 2p23pl2A 3 + 5-P33A3Q1333 

-P22P13A3 + 4pi2pi 3 + 2pi2pi 3 + 2pi 2 A 2 - P22A1 
+ 2P23P12A3 + 4p 33 Ai + 2p 33 pi2p22 - IP33A3Q1333 

Pl2 + 4p33A + 4p 33 pf 3 + f pl3<2l333 



2A 



2G 



14) 
14) 

14) 

15) 
15) 

15) 
15) 
16) 

16) 

16) 

16) 

17) 
17) 

17) 
18) 

18) 



Pl23 



Pll3p223 



Pll2p233 



Pl22pl23 



Pll2p223 



Pl22 



P123P222 = 2p22p23pl2 ~ 2pf 2 + 2pi 3 p| 2 - 2p 23 p3 3 Ai - 2A 3 pi 3 p 33 
+ 2p22pl2A3 + 8P33P13A2 — |pl3<5l333 + ^23^3*31333 

P122P223 = 4p22p23pl2 - 2pnp 22 + 2p^ 2 + 2p 22 p 12 A 3 + A§pi3p 33 
+ 4p23p33Al + P33A3A1 + 8p 33 A + |pl3<3l333 

— IP23A3Q1333 — 4p 33 p 13 A 2 — |A 2 Qi333 
P113P333 = 2p33p? 3 + 2 p| 3 pll ~ 2 Pl2 + 2p 3 3pl3A 2 - P23P33A1 

— 2p33Ao — |pl3Ql333 

P123P133 = 2p23p? 3 + 2 pl3Ai + 2pf 3 A 3 - 2p 23 A + 2p 33 pi 2 pi3 + \p\lQ 
P122P222 = 6p 33 p22Ai - P23A3A1 - 2p 33 pi 2 A§ - 2p 23 Ai + 2A3P13A2 

+ 2A 3 pnp 33 + p 23 pi 3 A 3 - 2pf 3 A 3 - |p22A 3 Qi333 

+ 4p 33 pi 2 A2 - 2A3A0 + 4p 22 pi2 - 3 pl2<3l333 - 8pi 3 Ai 

P113P233 = 2pi 3 A 3 + 2p 23 p? 3 + 2p 33 pi 2 A 2 + 2p i3 Ai + 2p 33 p 23 pn 

SP22A1 — 

Pll2p333 = 2p23pl3A2 - P23^1 _ 4pj 3 A 3 + 4p 33 pi2pl3 - 2p 23 A 



1333 



2p23Ao — P33P22A1 — |pl2Ql333 



+ 2p33p23pll - 2p23pf 3 + 3P12Q1333 

5P23P12A2 - |pllpl3 + ^Pl3p23pl2 - |p33pllp22 + IP33A2Q1333 
- 5P23P22A1 + IP22P13A2 + 3p 22 A + 5P12A1 + 5P23P11 
+ 3-P13P22 - Ip33^0 + |P33P?2 + 3pl2pl3A 3 ~ ip§ 3 A 3 Ai 



Pl22pl33 = -3PHP13 + f Pl3p23pl2 + f P33pllp22 



|pi 3 pn 



|P13P22 



|P33P?2 



5P23P22A1 - IP22P13A2 



+ 3pl2pl3A3 + IP33A3A1 + \p\2^-\ + 



3 

16,J 



P23P12A2 



3 P33^0 - |p33A2Ql333 
|pllPl3 + f Pl3p23pl2 + §P33P11P22 + f P23P1I + |p23pl2A 2 



+ |P? 3 P22 - §P33P? 2 



5P23P22A1 + |p22pl3A2 
3P33-^3Ai + 3P12A1 - |p| 3 A + |p33A2Ql333 



XPHPl3 + f Pl3p23pl2 + §P33P11P22 + f P23PII 



|p33A 2 < 



'1333 



|P? 3 P22 



|P33P? 2 



3-P23P22A1 - |p22pl3A2 - 6p22Ao 

+ f P23P12A2 + 3P33A3A1 + ^pi 2 Ai + ^p 33 A 
: 2p 33 A 3 A + 2pi2pi 3 p22 - 2pnpi2 + 2Aipi 3 p 33 + 2A 3 pf 2 
+ 2p2 3 p? 2 — 3-A1Q1333 — 3-P13A3Q1333 



P113P222 = -8P23P33A0 - 4p 33 A 3 A + 8Aipi 3 p 33 + 4P22P12A2 + |AiQi 333 

+ 4pi2pl 3 p22 + 4pnpi2 - A3P23P33A1 - 2p23pi2 + 2 p22p23pll 



4A 3 p? 2 - 2 p1 2 Ai + IP23A2Q1333 ~ f P13A3Q1333 



IP13A3Q1333 



8P33A3A0 + A 3 p 23 p 33 Ai - 4Aipi 3 p 33 + 2p 23 pi2 

+ 2p22p23pll + 8p23p33Ao + 2A 3 p? 2 - |p23A2<9l333 - 5A1Q1333 

-4A 2 p? 3 - 4A 2 A + Af - I6A0P13 - 8A pi 3 + 4p 22 p? 2 
^1 



+ 2A 3 pi 3 Ai + pi 3 A 3 + 4Aipi 3 p23 + 4Aipi2p 33 

+ 8p 33 p22Ao — |pl2A 3 Ql 333 + |pllQl 33 3 — 4A3P23A0 

P111P333 = -9A 3 p 23 A - 2pi 3 A 2 + A2P23A1 - 2A 2 pnp33 + 2pn<2i333 
- 10A p 23 + 6p 33 pi 3 pn - 10A pi 3 - 2pf 3 + 7Aipi 3 p 23 
+ 2Aipi2p 33 - 2p 33 p 22 A + 6A 3 pi 3 Ai + 2A? - 6A 2 A - 4A 2 pi 3 
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Pll3pl33 = 2p33pl3pll + 2A 2 pi 3 + 2 A pl3 + 2A p| 3 + 2pf 3 
+ Aipi2p33 - 2p33p22Ao - 2Aipi 3 p23 

P112P222 = 4A 2 pi 3 + 2P23A3P13A2 - 2A2P23A1 + 4A 2 pnp33 - P23A3A1 

+ 4A 2 A - 2A 2 + 8A pl 3 + 2pnp 2 2 + 2p 22 p 2 2 + IP12A3Q1333 

- 2A3P33P12A2 + 2A3P33P22A1 + 4pi 3 A| - 8A3P13A1 - lpllQl333 



Pll3pl23 — 



Plllp233 — 



8A1P13P23 + 4Aipi2p33 + 8A3P23A0 

„2 



2 \ 
3A2 



2p 2 ,A 2 



1333 - z Pl3 A 3 

-2P33A3A0 + 2A 2 pl2pl3 + 2pl2p? 3 + 2pi 2 A - 2p 23 p2 2 Ao 
+ 2pl 3 p23pll + 3-P33A1Q1333 

2pl2p33pll - 2pi2pi 3 + 6A1P12A3 + 6Aipi2p23 
+ 4pi3p23pll - 2pi 2 A 2 - 2pnAi - 4A 2 pi2Pl3 + A2P22A1 

- 2A 2 p23pn - 8P23P22A0 + P13P22A1 + 4p| 3 A 3 Ao 
+ 6piipi 3 A 3 - 9A3P22A0 + 2pi 2 Ao - IP33A1Q1333 
P112P133 = 4p 33 A 3 A + 2pi2p 33 pn + 2A 2 pi2pi3 - 2pi 3 p22Ai + 2pi 2 A 

+ 2pl2pf 3 + 4p 23 p22A - |p33AlQl333 - Aipi2p23 



P113P122 = 2p22pl3pll " 
+ 2pi2p22Al 

+ 8P13P33A0 + |AoQl333 + I P23A 



4p33A 2 A + 2p 2 2 pi 3 - 4p 2 . 2 A + P33A 2 

4p23p33A3Ao + A1A3P13P33 — IP13A2Q1333 



P112P123 — 



P111P223 — 



P33A2A0 + 2p 2 2 pi 3 + 2pi2p23pll + 2p| 2 A 



I A0Q1333 — 5P23A1Q1333 



P111P222 — 



P112P122 — 



P113 



1V1333 

2p 33 A 2 

+ 2A 2 p 2 2 ~~ 2pi2p22Al + 2p23p33A 3 Ao - 

-12p3 3 A 2 Ao + 2p 22 pi3pii + 6pnA 3 pi2 - 4p 2 j + 4p 33 A 2 

- 4pl 2 Ao - 2A 2 pnp22 - 8pi 3 p 33 A - f P23A1Q1333 

- A2p33A 3 Ai + 8p23p33A 3 Ao + |pl3A2Ql333 + §A|Ql333 

- 2A 2 p 2 2 - 2p 2 2 pi3 + 4pi2p23pll - AiA 3 pi 3 p 33 

- 2A 3 AiQi 333 + 4pi2p22Ai + 3A0Q1333 + 9p33A§Ao 

6A3A2A0 — 6A1A0 — 2A3A 2 — 2p22Ai<5i333 — 6A§pi 3 Ai — 8A3A0P13 

- 2pi2 + 2p 12 A 2 (5i333 + IOA3A0P23 + 4p 23 A^pi 3 - 2A3P11Q1333 
+ 6pnp22pl2 - 4p 33 pi 2 A| - 16A pi 3 p23 + 9A 2 p 2 3A - 4p 23 A 2 

- 2A 3 A 2 p 2 3 + 2A 3 p i3 A^ + 6A2P13A1 - 8A0P12P33 + 10Aip 2 3 
+ 2A 3 A 2 pnp33 - 2p 2 3 A 2 Ai + 4A 3 Aipi2p33 + 2A2P33P22A1 

- A 3 A2p23Ai + 2pnp 33 Ai - 7A 3 Aipi3p23 + 2A3A0P22P33 
2pnp22pi2 + 2pnp 33 Ai - 2AiA - 8A 3 A pi3 + 2pf 2 

- 2A 3 A p 23 + 8A pi2p33 + 2A 3 A p22p33 + 2A 2 pi 3 Ai 

+ A 3 Aipi 3 p23 - P23A? - 2Aip 2 3 - |pi 2 A 2 (5l333 

4p 2 3 pn - 4A pi3p22 - 4A pi2p23 - 4p§ 3 A 2 A + p| 3 A 2 



+ 4Aipi2pl3 + |p33Ao( 

P111P133 = 8A0P12P23 - A1P12A2 - 



/1333 

6A2P22A0 



- 2A 2 pnpi3 + 2p22A 2 
2p33-^i - 6Aopn + 2p 33 p 2 1 + 9A pi 2 A 3 - IP33A0Q1333 



'§ 3 A 2 A 



- Aipi2pl3 - 2Aipnp23 + 2p 2 3 pn - 4A pl3p22 - 
Pll2pll3 = -4p23p33A2Ao + 2Aip 2 2 + P23P33A 2 - 4A pl2p22 

4pl2pl3pll + f P23A0Q1333 - §Pl3AlQl333 + 4A A 3 pi3p33 
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(-23) 



(-24) 



(-24) 



(-26) 



(-27) 



(-30) 



PlllPl23 



PlllPl22 



Pll2 



PlllPll3 



Plllpll2 



Pill 



3 

' 3P23A0Q1333 + 3P13A1Q1333 

-3A3A1A0 — A2AI + 4A^Aq — 8Ag 



?1333 



2pl2pl3pll - 2p33A 3 Af + 2Aip 33 A - 2A A 3 pi3p33 

+ 2piipi 2 A 2 - 2pnp22Ai + 6A2P33A3A0 + 8p2 3 p 33 A2A 
+ 2A0P12P22 - 2P23P33A? + 3-A2A1Q1333 - 3A 3 A <9 

2p23p?l 
AlP? 3 ^3 

+ 4A2A0P23 — f P22A0Q1333 + IP12A1Q1333 - 



8A1A0P23 - 2p 2 3 A 2 



§A2pllQl333 + Af P13 



+ 8A 2 A pi 3 - 9A§A pi3 - 2A 3 A 2 P23 + 2pf 1 p 22 - 2A 2 Aipi 2 p33 

+ 2pnpi2 + 8A0P13 + A2A3P13A! + 4A 3 A pi2p33 + 2A 2 Aipi 3 p23 

- 8A3A0P13P23 + 4A2P33P22A0 + 2A 3 Aipnp33 + 6A2A3P23A0 
: -8A0P13 + 8A pnp33 - 4AiA p23 - 4A 2 A p2 3 - 8A 2 A pl3 

+ 4 PHPl2 - 8A + 4A 3 A pl3p23 + 4A 3 A pl2P33 + p| 3 Af 
+ 4A^pi3 + |p22AoQl333 — IP12A1Q1333 

: 4A0P33P13A2 + 4pi!pi 3 - P33A2A? + 4p 33 A 2 A - Af pi3p 33 

- 3p33hh\) + 2pi2pnAi - 6p2 2 A pii + 2A p 2 2 



2A2A0Q1333 — |Aopl 3 Ql333 + §AiQl333 + 8p33Ap 



16p 23 Ag 



: — f P12A0Q1333 + 2A3pi 3 Af + 8A2A1A0 + Af P13P23 — fpnAi 
+ 6A3P33A0P11 + 4p 2 1 p i2 + 4p 23 A^A - 18A 3 A^ - 2Af 
_ 2A 3 p^ 3 Ao - 6A3P13A2A0 - 3A1A3P23A0 + 2A1P33P22A0 

- A2P23A? - 2Aipi2p33 - 2A1A0P23 + 8A1A0P13 + 4A p33pi2A 2 

: -4p 23 A? - 8p^ 3 A^ - 8pi 3 A|A - 4A 2 A p 2 3 + 16p2 3 A 2 AiAo 

+ 4p23AiA pl3 + 6A3A1A0P13 - 4A pnQl333 - 36A3P23A0 

+ I8A3A2A1A0 + A|Af - 4A 3 A? - 4A|A - 4A 2 A - 27A§A§ 
+ 12A 2 A^ + 2pi 3 A 2 A 2 + 16pi 3 A§ + A 2 p 2 3 - 4pi 2 Aip 3 3A 



1333 



4pllp33Af + I6P11P33A2A0 + 4pf x + 8p 3 3p2 2 A^ 



E The three-term three- variable addition formula 

The addition formula in Theorem 15.31 is constructed using the following polynomials. 



- jsd 3 Q 1333 (u)p^(v)d 3 Q 1333 (w) + ^d 1 Q 1333 (u)p^(w)p 133 {v) 

- 3 |pll2(M)pll2('y)p [111 (^) - j 2 -p22(v)d 1 Q 1333 (ll)d 1 Q 1333 {w) 

+ M 5 3Ql333Np [111 (-")pm(™) + ^P11 2 (M)P [221 (-")P22 2 (^) + § p [23] («) p 2 2 («) p [23] (w) 
+ lP333(M)plll(^)p [22l ('y) + fQl333(M)p [33l («)<3l333(™) + |pn 2 («)pl3 («)pll2 (w) 

+ ^d 2 Qi 333 (u)d 2 Qi 333 (v)p^(w) - |p222(M)piii('w)p [12l ( u ') 

+ ^d 1 Q 1333 ( U )p 12 (v)d 2 Q 1333 (w) + ^(u)p^( V )p^(w) 

+ |plll(M)pl22(^)p [11I (' U ) - |P [331 ( M )P123(V)P12 3 ('W) + | P [33] («) P223 (v) Pll3 (w ) 

- lQi 333 (u)p 113 {v)p 113 {w) - ^Pi2 2 (m)p [221 (^)Pi2 2 (^) - %p [23] (u)p 12 (v)p^(w) 

- jdiQ 1333 (u)p u (v)p 133 (w) + pW(u)pW(w)p 13 {v) + lp^{u)p ll {w)p^(y) 
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§p[ 33 ](M)p 33 (>)p [33l (>) - Jsd2Ql333(u)d 2 Ql333(w)pn{v) 

jP ll2] ( u )p ll2] ( v )p l22] ( w ) - ^p in (v)d 3 Q 1333 (w)p 13 (u) + ipii(tt)pii(w)p u (w) 
ip[ 23 ](M)a 2 Q 1 333(v)p 1 3 3 (w;) + \piu(u)p 122 (w)p 13 (v) - |p 23 (m)Piii (w)pii2 (v) 
hp233(u)p l23 \v)diQi 333 (w) + ^p^(u)pn(v)p 12 (w) - lpin(u)pui(v) 



P 27 (u,v,w) = ^Pi2(v)9iQi 333 (m)Pi 33 (^) + lpn(w)p 133 {v)d 2 Q 1333 (u) 

- TeP ll3] ( w )pi33(v)d 2 Q 1333 (u) + %pn 3 (w)p 123 (v)Q 1333 {u) - §pl 22 l 0)p [22] (v)p 23 (u) 

+ ^ d 3Ql333(w)d 3 Q 13 33(v)p^(u) - ^p^\v)p 112 (w)p 333 {u) - ^ (v)p^ (w)p 13 (u) 
+ ^P23(w)p 112 (v)p 112 (u) + f p[ 12 ](™)p [12l ("u)p [12l (M) + ^Pl 22 (w)p 122 (v)p^(u) 

- j E pi 22 (v)pi 3 (w)p 1 i 2 (u) + j E d 3 Q 1333 {v)p^(w)p 122 (u) - lp 12 (u)p^\w)p 12 {v) 



Iplll('w)p333(^)p [12l ('« 



pii2(v)d 3 Q 1333 {w)p 13 (u) 



P 24 (u,v,w) = \p 13 {u)pW{v)pW{x 



8p222(«)p [12] (u)9 3 Qi 333 ( 



^lOi 333 (M)p 12 (w)p2 33 (l«) - Apl22(M)pll2(w)p23(w) + |p [12] Hp 333 (^)pll 2 ('w) 



^[ 12 ](«)p23(«)p[ 22 ](^) 



8' 

lp 12 {u)p 12 {v)p^\w) + ^p 222 {u)p 112 {w)p 13 {v) 
\d 3 Qi 333 {u)pi 22 (y)pi 3 (w) - fpi 33 (u)pi 33 (t!)p[ 13 ](i(;) + ^Pi 22 (m)pi22(v)pi3(^) 

^9 3 Qi33 3 (M)pii2(f)p23(^) + JPl33{u)p l33 (w)pn{v) - jQl 333 (u)p 123 (v)p 123 (w) 
j g p lll] (v)d 3 Q 1333 (w)p 122 (u) + !pll 3 (u)p33(™)pll 3 (f) + f Pl2 (v)p [23] (w)p 22 (u) 
kpl33(u)pl 2 (v)d 2 Q 1333 (w) - 1 ^P122(M)P [121 (' U )P222(^) + \ Pl3 (« )p [12] O) p [12] (w) 

jp33(u)p^(w)Q 1333 (v) - lp 223 (u)p 113 (w)Q 1333 (v) - -^pi 22 {u)p 122 (v)p^(w) 

3 
4 



Pi2Hpi2(w)pn(w) + 3p 333 (u)pm(>)pi 3 (™) - |p [llI ( ^ ')p [12] ( ^t )p [12l ( ^t, ) 

^lPl3(u)d 3 Qi 333 {v)d 3 Q 1333 {w) - 5§8 p'" 1 ( M ) 5 3Ql333('w)9 3 (5l 3 3 3 (w) 



2l^3Ql333(M)p333(«)p [22l ( w ) ~ fplsHplsMp' 221 (™) + \p 333 {u)p^ (v)p 122 {w) 



\pW(u)pW{v) + ±p 112 {u)p ll2 {v) 



A 2 + 



- 4pll3(«)p33(«)pll3(™) 



4pi3 3 (M)pi3 3 (lu)pii(w) + ■^Pi33(w)pi3 3 (w)p[ 13 ](w;) 



|pi3(m)Pi3(^)p [221 (w) 



\2 



P 21 (u,v,w) = %p 233 {u)p^(w)p 133 (v) - lpn 3 {u)p 123 (w)p 33 (v) + lpi 2 {u)p l2 {v)p 12 (w) 
+ l±pM(u)p 23 (v)p 13 (w) - \pM(u)p 13 (w)pM(v) - \pW(u)p 23 (v)pW{w) 

- lp [2?,] {u)p 22 {v)p 22 (w) + ^pl 22 (u)p 122 (w)p 23 (v) - |p[ 13 ](M)p 22 (v)pl2(^) 
+ Ig93Qi333(M)p[ 11 l(li;)p222(v) + ^<93<9l333(w)p23(»pl22(u>) 

§ed 3 Q 1333 (u)p 23 (w)d 3 Q 1333 (v) - ^pll 2 ('w)pl22(^) - ^P 1 " 1 («)pll2 (™)p333 0) 



|P22H9 2 Q1 333 (W)P1 3 3W + ip[ n ](«)p [111 (^)p [121 ^) 

\p 11 (u)p 233 (w)p 133 (v) - gp333(w)plll(t>)p23(u>) ~ jQpll2(u)d 3 Q 1333 (v) 
|P23(M)P [121 (^)P [121 (^) + TOp233(tt)9lQl333(^)p22(v) 

|P333(m)P122(«)P [121 (^) + ^p233(M)9 2 Ql333(«')pl2('y) + f P [22] (v)p^ (w) 

lpn 2 {v)p 13 {w)p 333 (u) - j E p 222 {u)p 13 {w)p 122 {v) + lp^ 2 \u)pi 3 {v)p 13 {w) 
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|P22(«)pll(«)pl2(«>) + TeP222{u)p [11] (w)p 122 {v) - §pi 2 (li)pl33 (v)p 133 (w) 
J E p [22] (u)p 333 (v)p 22 2(w) + J S p222{u)p 13 (w)d 3 Q 1333 {v) 

!:Ql333(u)pl2 3 (v)p22 3 {w) + ± p 2 22 («) P222 (w) p [12] (v) 



Ai 



P [22] («)P23(«)P13( 



pl 3 (u)pl 3 (w)p [12] (™) + jW 12l C")pl22(w)p333(l>) 



JaPll2{v)p 13 {u)p 333 {w) + ^<93Ql333(M)p333(«)p [121 («?) 



2^3 



P 18 (u,v,w) 



- lpW(v)pW(w) - lpW(u)pW(v) + ^p ll2 {v)p 222 {u) 

- §Pi3{u)p2 3 {v)p [12] {w) + p 333 {u)p 13 (v)p 122 (w) + 2p 13 (u)p 13 (v)p 13 (w) 

- '%p22(u)pl 2 (v)p 12 (w) + lpl 3 (u)d 3 Q 1333 (v)p 333 (w) + ^pl 22 («)pl 22 (» 

- T{-d 3 Q 1333 (u)d 3 Q 1333 {w) + lp 233 (u)p 133 (w)p 12 (v) - \pt- 22 Xu)p 23 {v)p 23 (w) 

- Ip2 33 {u)p 233 {w)p^{v) + \p^ 22 \u)p 333 (v)p 333 {w) + |pll 2 (w)p23(w)p333W 

- |P33(m)P [33] (' u )P33(^) + fp233(w)p233(«)pll(™) ~ f p223(")p223(»<2l333(K') 

- ^p [111 (M)p22 2 (u)p222(™) + f P13(M)P222(«)P22 2 ('W) + 3pi 23 («) Pl23 («) P33 {w) 

+ lp ll3] (u)p 22 (v)p 22 (w) - ^p 222 (u)p 23 (v)d 3 Q 1333 (w) - lpn{u)p 22 (v)p 22 {w) 
-3p[ 11 ](«)p 13 («)pi 3 (t«;) + ipi 22 Np333Wp [111 M-iP [111 («)p [111 WP [111 W 

\p [11] (u)p 333 {v)d 3 Qi 333 (w) + §p 223 (»Pll3(«)p33(™) 
P233(m)P22(™)9 2 <2i333('w) + kQl333{u)Ql333(w)p 33 (v) ~ Jplll(u)p 333 (w) 



+• %Pi33{u)pi 33 (v)p 22 (w) + jp [12] {u)p 333 (v)p 222 {w) + lp^{u)p 23 {w)p^ 12 \v) A 

3 iP2 2 (M)pl 2 (^)pl 2 (^) + pl 3 (u)p 13 (v)p 13 (w) - ^pl 22 (M)<9 3 Ql333(lf) 
?1333(M)Ql333('y)p33(^) + jpl 3 (u)p 23 (v)p^ (w) + ±pl 12 l (w)p [121 (w) 
+ Ip333{u)pl3{v)pl 22 (w) + ±p 13 (u)p^(w)p 13 (v) + ^pl 3 (u)d 3 Q 1333 (v)p 333 (w) 

- ^P l22] (u)p 13 (v) - -^d 3 Qi 333 (u)d 3 Q 1333 (w) - ygpi 22 Gu)pi 22 (™ 

Jp [22] {u)p 13 {v) + ^P122(m)P122(«) - j3gP33 3 (u)p[ 22 ] (v)p 333 (w) 

jpi3{u)p 23 {v)p [12] {w) + ±p 13 (u)p 333 (w)p 122 (v) + \d 3 Qi 333 {u)pi 22 (w) 

Jg,d 3 Qi 333 (v)d 3 Qi 333 {w) - ^<9 3 (2i333(w)p333(u>)pi3(l>) - §p [12] (u)p [12I ( 
|pl23(v)pl 2 3(w)p33(M) + ^p333(w)p [12l Mp22 2 (t') 



Ai 



- \pl3{u)p 13 (v)p 13 (w) + ^pU 3 (w)p 33 ('U)p223(«') + f p23 (u)p [22] {w)p 23 (v) 

+ Tep23{u)pii 2 (v)p 333 {w) - ^p 33 (u)Q 1333 (v)Q 1333 (w) 

+ JeP233(u)pl33(w)pi 2 (v) - |pl 2 (w)pl2(«)p22(w) + § P33 (u)p^ (w)p 33 (v) 
+ ll;p333('")pl 22 (w)p [11I («') + j8d 3 Q 1333 (u)p 333 (w)p [11] {v) A3A1 



P 15 (u,v,w) 



Pl3{u)p ll2] (w) - ■^p222{u)pl 22 (v) - %pl 33 (u)p 22 (w)p 233 (v) 



|p[ 12 ](M)p 23 (t>)p 23 (u;) - ip 3 33(M)p33 3 (v)p[ 12 l(u;) + ^p 333 (u)p 112 (v) 
!P [121 (V)P [111 (™) - %p 333 (u)p 23 (w)pi 22 (v) + ^p 22 {u)p 22 {v)p 12 (w) 

Tsp222(u)d 3 Q 1333 (w) + lp 13 (u)p 13 (v)p 23 (w) - ^p 333 {u)p 13 (w)p 222 (v) 

f P233(m)P233(«)P12(» ~ f P33(«)p223(u>)pl23(w) 
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^rp222(u)p [11] (v)p 333 (w) + f p 2 z{u)p [11] {w)p 13 {v) 



A3A0 



JPl3{u)p ll2] (w) - \p^(u)p 13 (w)p 23 {v) - ±pl22{u)p 23 {w)p 333 (v) 

+ JP l22] (u)p 23 (v) - jgp23{v)p 333 (w)d 3 Qi 333 (u) + ygp222Npl22(» 

- Ip23(u)p23(v)p ll2] {w) - \p 13 {u)p 13 {v)p 23 {w) - Jq p 1 3 (v)p 222 (w)p 333 (u) 



JP [12] (u)p^(v) + j g p222(v)8 3 Qi 333 (w) - §p22(»p22(v)pl2(» 



A2A1 



Pi3(u)pi 3 (v)p 23 (w) + i2p 3 33(u)p33 3 (u;)p[ 12 ](t;) 



A3A1 



P 12 {u,v,w) = |pi 3 (t;)pi3(u;) - ^p^(u)p^(v) - jp 333 (u)p 13 (w)p 333 (v) 

+ Jgp222(u)p222(v) + § p23 («) pl 3 (to) p 23 (v) - § p 2 2 (u) p 2 2 (v ) p22 (w) 

+ ^P333{u)p lll] (w)p 333 (v) + \p^{w)p l3 (u) + ^d 3 Q 1333 (v)p 333 (w) 

+ fpl22(«)p333(» ~ §P23(m)P [111 (' u )P23(^) + f P233 0)p22(™)p233 0w) 
- JP23{U)P [12] ( W ) + |P333(")P222('W)P23(W) + § P223 (u) p 3 3 (w) p 2 23 («) 



4p 33 (it) p33 («)<5l333 O) 



A2A0 



p33(M)p33(^)<5l333(l') 



^Pl22(u)p333(H 



49 3 0i333(v)p333('U;) - -fp 23 {u)p 13 {w)p 23 {v) 



16 

M P333 (u) P [11] (w)p 333 (v) 



\p333(u)pl 3 (w)p 333 {v) 



\ 3 A 



J2p33(v)p223{u)p223(w) ~ ± p 2 22 («) P222 («?) ~ ^ Pl3 (») Pl3 («>) 



P23(")pl3(^)p23(^) - ^P333(m)P [11] («)P333(^) + | P333 («) Pl3 («) P333 (w) 
-Pl22(u)p 3 33(^) + ip [22] (^) - 4P22(m)P233 0)P233O) 



+ f P33(«)p33(™)Ql333(«) - ^9 3 Qi33 3 (t))p333(u;) + §p 2 2(u)p22(w)p2 2 (™) 



\p2 3 {v)p^{u)p 23 {w) + \p^(v)pW(w) - \p^{v)p 13 {w) 
\p23{u)p222{v)p 333 {w) - ip2 3 Gu)p [12l O) X{ + f Pl3 C«) Pl3 ( 



to 



+ JsP333(u)d 3 Q 1333 (v) - ^p 33 (u)p 33 (w)Q 1333 (v) + Ygpl2 2 K)p333 0) 

+ |p23(M)pl 3 (w)p23(«) - 1^P333(W)P13(^)P333('W) A 3 A 2 Ai - \p 13 (v)p 13 (w)Xl 



P 9 (U,V,W) = 32p23(M)p333(«)p333(^)A 3 Af - |p 23 («)p23 («)p23 («>) A3 A 2 

- |pl3(«)p23(^)A 3 A 2 Ao - |p23(M)p23(u)p23(l«)A 1 Ao - ±p 2 22 («)p333 (™) A3 A 2 

- |p23(w)p333(v)p333(^)AlAo + §p333(w)p222(v)A 3 A2Ao - fpl 12 ! (t?)Ai A 

+ \p [11] {v)p 23 {w)\ l \ + lp23(u)p23(v)p23(w)\ 3 \ 2 \o + f Pl3 (™)p23 («) A3 A? 
+ 3Pl3(w)p23(«)A|Ai + xPl3(™)p23(«)AlA + ^p222(lt)p333(«)AlA 



P 6 (tt,t;,t«) = 9p 13 {w)X% + %p 33 {u)p 33 (v)p 33 (w)XlXl + ^p3330)p333(™)A§A 2 Ao 



P [111 ( M ) A + §P33(M)p33('w)p33(lf)A2Ao + ^p333(^)p333(™)A 3 Ai A 
P23(m)P23(«)A 3 AiA - 27p 23 (t>)p2 3 (^)Ao - §pi3 Cu>)AiA 
P333(W)P333(«)A1A + 7;pl 3 (™)A 2 A 2 - |p333('")p333(^)Ao 



+ 6p 33 (M)p33(»p33(™)A 2 | - §P33(m)P33(«)P33('w)A2A 2 + |p23(v)p23(w)A|Ao 
~ !P33NP33 0)P33(™)A 2 A 2 A - 3p 33 {u)p 33 (v)p 33 (w)X 3 X 1 X 
- JeP333(u)p 333 (v)XlXl - |p23(M)p23(«)A2A? + ^p333O)p3330w)A 2 A 2 
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P 3 (U,V,W) = -±±P2 3 (W)\ 3 1 - f-p23(«)A 3 A§ + ^p23(w)\ 2 X 1 X 

Po(u, v, w) = IA3A2A1A0 + gA|Af — ^3X1 — 5A2A0 — |Af Ao — ^-A|Ag + |A2Aq 
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